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ABSTRACT 

Propagation  of  electromagnetic  waves  through  a  homogeneous  aniso¬ 
tropic  column  of  a  medium  of  Infinite  length  Is  considered.  The  anisotropy  of 
the  medium  Is  characterized  by  the  dyadic  form  of  the  permittivity  4  and  per¬ 
meability  ju  .  This  anisotropic  column  Is  surrounded  by  a  coaxial  homogeneous 
Isotropic  medium  characterized  by  scalars  £  3  and  ^  ,  this  complete  structure 
being  enclosed  by  a  perfectly  conducting  metallic  circular  cylindrical  wave¬ 
guide.  A  magnetic  current  ring  source  Is  Inserted  symmetrically  In  the  Iso¬ 
tropic  medium.  Throughout  the  analysis  the  strength  of  the  source  Is  con¬ 
sidered  to  be  an  arbitrary  function  of  the  polar  angle  0.  For  this  general 
problem  the  complete  expressions  for  fields  (due  to  the  source),  power  flow, 
and  the  dispersion  relation  have  been  studied. 

To  solve  a  source  problem,  dyadic  Green's  functions  for  both  point 
electric  current  source  and  point  magnetic  current  source  have  been  con¬ 
structed  In  a  formal  way  from  the  source  free  solutions  of  the  appropriate 
Maxwell's  equations.  These  dyadic  Green's  functions  can  be  used  for  any 
arbitrary  source. 

From  the  general  expressions  for  the  transverse  fields  In  terms  of  the 
longitudinal  fields.  In  any  arbitrary  cylindrical  region  (unbounded)  the  propaga 
tlon  wave  number  of  a  TEM  mode  travelling  In  the  longitudinal  direction  z  has 
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been  obtained.  The  wave  numbers  for  a  TEM  wave  propagating  In  a  direction 
perpendicular  to  z,  can  be  obtained  from  the  general  expressions  for  the  trans¬ 
verse  wave  numbers,  using 


b_ 

bz 


-jX  =  o  . 


The  results  of  the  above  general  problem  have  been  used  to  study  the 
wave  propagation  In  an  anisotropic  plasma  column  and  an  anisotropic  ferrite 
column  separately.  The  various  possible  pass  bands  for  the  propagation  of 
electromagnetic  waves  In  an  anisotropic  plasma  column  have  been  obtained  and 
a  special  case  Is  considered  for  numerical  computation  of  the  longitudinal 
electric  field  In  a  plasma.  The  analysis  for  the  plasma  problem  emphasizes  the 
slow  wave  propagation. 
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I 

INTRODUCTION 

Electromagnetic  wave  propagation  through  an  anisotropic  medium  has 
been  studied  by  many  authors  In  different  situations  -  from  Ionosphere  to  con¬ 
ventional  waveguides.  Ionospheric  anisotropy  Is  due  to  the  presence  of  the 
earth's  static  magnetic  field.  Ionized,  but  macroscoplcally  neutral  gases  of  any 
kind  are  known  as  plasmas.  An  Ionized  and  neutral  stationary  plasma  In  a  weak 
electromagnetic  field  can  be  represented  as  an  equivalent  dielectric  medium. 
Moreover,  when  this  plasma  Is  situated  In  a  uniform  static  magnetic  field,  the 
strength  of  which  Is  not  necessarily  small.  Its  equivalent  dielectric  "constant" 
behaves  as  a  dyadic  (tensor).  It  Is  well  known  that  radio  wave  propagation 
through  the  Ionosphere  depends  on  the  frequency  of  the  electromagnetic  wave, 
electron  density,  the  collision  frequency,  the  ton  gyrofrequency  and  the  electron 
gyrofrequency.  The  same  Is  true  for  a  wave  propagating  through  a  plasma 
waveguide.  This  knowledge  of  the  wave  propagation  Is  essential  for  sattsfactoiy 
long-distance  radio  communication  through  the  Ionosphere. 

The  study  of  a  plasma  Is  also  of  vital  Importance  In  the  field  of  thermo¬ 
nuclear  reactions.  In  a  thermonuclear  reactor  a  static  magnetic  field  Is  used 
to  confine  and  to  heat  a  plasma.  Many  times  It  Is  desirable  to  obtain  Information 
on  the  temperature,  density,  etc.  of  the  plasma.  Such  Investigations  which 
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determine  the  characteristics  and  behavior  of  the  plasma  are  known  as  "plasma 
diagnostics. " 

Besides  the  above,  the  propagation  of  electromagnetic  waves  In  plasma- 
filled  or  partially  plasma-filled  waveguides  has  aroused  considerable  Interest 
In  recent  years  ^6^  [  7  J  ^8^  ,  primarily  because  of  possible  applications  to 
the  generation  or  amplification  of  microwaves.  A  medium  whose  dielectric 
"constant"  Is  a  tensor  Is  called  gyroelectrlc. 

On  the  other  band,  there  Is  another  class  of  materials  known  as  ferrites 
which  exhibit  ferro -magnetic  properties.  The  chemical  composition  of  the 
ferrites  may  be  expressed  ^3  J  £is]  by  the  formula  MO Fe^ Og ,  where  M 
represents  a  metal,  such  as  Mn,  Fe,  Nl,  Cu,  Mg,  Al,  Co,  etc.  Although 
ordinary  Iron  (Fe)  and  nickel  (Nl)  possess  ferromagnetic  properties,  they  are 
of  little  use  as  microwave  components  due  to  their  high  losses.  But  the  ferrite 
materials  mentioned  above,  whose  specific  resistances  are  above  10®  times 
higher  than  those  of  the  metals,  with  relative  permeabilities  ranging  up  to 
several  thousands  and  rel.  we  dielectric  constants  varying  from  5  to  25,  have 
extensive  use  In  microwave  devices.  In  the  presence  of  a  static  magnetic  field 
the  permeability  /u  of  a  ferrite  becomes  anisotropic,  l.  e. ,  /u  becomes  a  dyadic, 
which  Is  the  characteristic  of  a  gyromagnetlc  medium.  The  medium  whose 
dielectric  constant  and  permeability  both  are  tensors.  Is  known  as  gyrotroplc. 
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Id  connection  with  the  Faraday  rotation  of  guided  electromagnetic  waves 

In  a  gyromagnettc  medium  with  a  uniform  static  magnetic  field  In  a  circular 

cylindrical  waveguide  (the  axis  of  the  waveguide  coincides  with  the  gyro-axis 

which  Is  also  the  direction  of  the  static  magnetic  field  and  has  been  taken  as  the 

zaxls),  Suhl  and  Walker  j^2  J  have  shown  that  only  circularly  polarized  modes 

exist.  If  Ez  j  0,  and  Hz  j  0,  and  pure  TE  and  TM  modes  do  not  exist. 

However,  pure  TE  and  TM  modes  can  exist  If  - —  =  0.  It  should  be  noted 

dz 

that  In  this  case  as  well  as  throughout  the  present  work,  the  anisotropic  medium 
under  consideration  Is  homogeneous.  For  an  anisotropic  plasma  medium  TE 
and  TM  modes  can  exist  £  7J  Independently  In  another  special  case  when  the 
axlal-statlc  magnetic  field  Is  Infinite.  Besides  the  above  mentioned  work,  a 
number  of  Investigators  Including  Van  Trier  £  1 J  ,  Gamo  j^ldj  ,  Falnberg  and 
Gorbatenko  [4J,  Agdur  ,  Epstein  £3]  ,  Trlvelplece  ^8j  etc.  have 
carried  out  research  In  connection  with  wave  propagation  In  gy roe lec trie,  gyro- 
magnetic,  or  gyrotroplc  media  with  various  configurations.  All  of  the  research 
work  cited  above  except  that  In  (^7  J  considered  only  the  source  free  resonance 
behavior  of  electromagnetic  waves.  In  die  present  problem,  however,  a  source 
of  electromagnetic  waves  which  Interact  with  the  anisotropic  medium  Is  Included. 
In  Appendix  A,  a  general  formulation  of  the  source -free  problem  Is  presented 
for  any  cylindrical  geometry  with  arbitrary  cross  section.  This  formulation  Is 
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suitable  even  for  an  unbounded  anisotropic  medium,  provided  cylindrical  sym¬ 
metry  Is  assumed.  In  Appendix  B,  dyadic  Green's  functions  for  a  point  souroe 
(electric  current  or  magnetic  current  or  both)  are  constructed  from  the  general 
source  free  solutions  of  the  Maxwell's  equations  for  both  dissipative  and  non- 
dissipative  media.  For  such  a  construction  of  Dyadic  Green's  functions  refer¬ 
ences  ,  £llj  ,  and  ^12j  have  been  found  very  useful.  An  alternative 
method  using  a  transmission  line  formulation  can  be  devised  for  the  construction 
of  Dyadic  Green's  functions. 

In  chapter  I  the  problem  considered  Is  to  find  the  dispersion  relations  and 
the  complete  fields  due  to  an  excitation  by  a  magnetic  current  ring  source  situated 
In  a  cylindrical  Isotropic  homogeneous  medium  characterized  by  a  relative  di¬ 
electric  constant  £  a  and  a  relative  permeability  ^3  ,  which  encloses  a  central 
cylindrical  column  of  a  homogeneous  anisotropic  medium  characterized  by 
dyadlcs  C  and  £  ,  this  whole  structure  being  enclosed  by  a  perfectly  conduct¬ 
ing  cylindrical  waveguide.  This  general  analysis  has  avoided  specifying  any 
particular  medium,  say  a  plasma  or  a  ferrite,  and  also  It  does  not  necessarily 
consider  a  ring  source  of  constant  strength.  A  ring  source  (magnetic  current) 
represents  an  Idealization  of  a  possible  excitation,  for  example,  a  circumferen¬ 
tial  slot  In  the  waveguide  wall,  or  an  annular  slot  on  a  thin  metallic  disc*  fitting 

*  In  this  example  one  must  also  consider  the  boundary  condition  for  the  conducting 
metallic  disc. 
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tightly  across  the  waveguide.  Although  a  ring  source  Is  taken  for  analysis,  any 
other  type  of  source  can  also  be  handled  adequately  since  the  formal  expressions 
for  dyadic  Green's  functions  for  an  electric  current  source  and  a  point  magnetic 
current  source  are  given  In  Appendix  B.  A  magnetic  current  ring  source  Is  more 
appropriate  for  a  plasma  problem,  whereas  for  a  ferrite  problem  an  electric 
dipole  at  the  center  of  a  cross  section  of  a  circular  waveguide  Is  more  appro¬ 
priate. 

In  Appendix  C,  the  general  dispersion  relation  of  Chapter  I  has  been 
evaluated  In  a  number  of  Interesting  special  situations  with  appropriate  limiting 
processes.  Although  In  these  special  cases  the  procedures  are  also  applicable 
to  obtain  expressions  for  the  total  fields  due  to  the  source  from  the  general  ex¬ 
pressions  given  In  Chapter  I,  no  attempt  has  been  made  to  obtain  these  expres¬ 
sions  owing  to  the  laborious  task  they  Involve. 

Chapter  II  deals  with  a  problem  In  which  the  anisotropic  column  Is  taken 
to  be  a  plasma  In  an  axial  static  magnetic  field,  using  the  results  of  Chapter  I. 

In  this  case  necessary  conditions  for  slow  wave  propagation  (which  give  maxi¬ 
mum  passbands)  have  been  obtained.  The  sufficient  condition  and  hence  the 
actual  passbands  can  be  obtained  from  the  solution  of  the  dispersion  relation. 
Since  It  Is  not  possible  to  study  a  dispersion  relation  In  general,  a  few  special 
cases  have  been  discussed.  In  addition  to  these  a  more  general  dispersion  re  la- 
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tlon  for  slow  wave  propagation  has  been  considered  for  numerical  computation 
In  Appendix  D.  In  this  case  the  lowest  eigenvalues  are  found  and  the  corres¬ 
ponding  longitudinal  electric  field  Is  calculated. 

hi  Chapter  in,  the  results  of  Chapter  I  have  been  applied  to  study  wave 
propagation  through  a  ferrite  column  with  a  uniform  axial  static  magnetic  field. 

It  may  be  noted  here  that  results  for  a  ferrite  problem  can  be  obtained  by  using 
duality  on  the  corresponding  results  for  a  plasma  problem  when  the  boundary 
conditions  on  H  In  the  ferrite  are  the  same  as  those  on  E  In  the  plasma,  for 
example  an  unbounded  plasma  and  a  ferrite. 

From  the  formal  expressions  of  the  transverse  electric  and  magnetic 
fields  as  functions  of  the  longitudinal  fields  Ez  and  Hz  ,  and  using  the  expres¬ 
sions  for  the  transverse  propagation  wave  numbers  obtained  In  Appendix  A, 
conditions  for  TEM  wave  propagation  In  the  direction  parallel  to  or  perpendi¬ 
cular  to  the  d.  c.  magnetic  field  have  been  obtained.  These  conditions  provide 
expressions  for  the  propagation  wave  numbers  In  the  respective  cases.  The 
conditions  which  give  the  possibility  of  a  TEM  wave  propagation  In  an  unbounded 
medium  cannot  be  valid  In  a  bounded  medium  or  In  waveguides  (except  those 
bounded  by  two  non-connecting  metallic  boundaries).  Therefore  these  discus¬ 
sions  suggest  that  the  study  of  a  dispersion  relation  for  a  bounded  anisotropic 
medium,  under  the  condition  of  TEM  wave  propagation,  Is  meaningless  and 
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Inconsistent.  It  can  be  shown  that  the  condition  of  TEM  wave  propagation  In 
the  direction  of  the  static  magnetic  field.  Is  equivalent  to  zero-value  of  the 
product  of  the  two  transverse  wave  numbers.  However,  results  for  such  a 
situation  have  been  presented  In  the  literature  mistakenly.*  A  reason  for  such 
Inconsistent  results  may  be  that  the  authors  overlooked  the  direct  equivalent 
relations  between  the  conditions  of  TEM  waves  travelling  In  the  direction  of 
the  static  magnetic  field  and  the  vanishing  conditions  of  the  product  of  the  trans¬ 
verse  wave  numbers. 


*Agdur,  6  pages  183  to  185,  also  ^4^ 


page  497. 
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I 

GENERAL  PROBLEM 

Statement  of  the  Problem 

An  Infinitely  long  column  of  an  anisotropic  medium  characterized  hy 
tensors  6  and  n  ,  of  radius  a,  Is  situated  coaxially  Inside  a  perfectly  con¬ 
ducting  circular  cylindrical  waveguide  of  radius  b.  The  annular  space  between 
the  column  of  the  anisotropic  medium  and  the  cylindrical  waveguide  Is  filled 
with  an  Isotropic  medium  characterized  by  scalars  and  ^  •  The  electro¬ 
magnetic  fields  are  Introduced  Into  this  system  by  a  magnetic  current  ring 
source  of  radius  c,  the  center  of  which  lies  on  the  axis  of  the  waveguide,  such 
that  a  ^  c  4  b.  The  total  fields  and  their  behaviors  are  studied.  Figure  1  shows 
the  geometry  of  the  problem. 

General  Formulation  of  the  Problem 

For  convenience  the  plane  of  the  ring  source  will  be  chosen  as  z  =  0, 
where  the  axis  of  the  cylinder  lies  along  the  z-axls.  Due  to  the  cylindrical 
symmetry  of  the  structure,  cylindrical  coordinates  r,  9,  and  z  will  be  used 
here.  If  the  ring  source  Is  very  thin  both  In  the  radial  and  In  the  axial  direction. 
It  can  be  represented  In  the  following  way 

I  =  9  I  e^  =  9  m(9)  S(r  -  c)  S(z)  e^ut  (1) 

— m  — o  in  -o 
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where  u  =  exciting  angular  frequency 

6  =  unit  vector  In  the  ©-direction 

-o 

m(d)  =  strength  of  the  source  In  volts 
Sir  -  c)  and  S(z)  are  well  known  Dlrac-delta  functions. 


FIGURE  1 

Although  ultimately  m(d)  will  be  chosen  as  a  constant  for  numerical 
computational  facility,  the  present  formulation  of  the  problem  Is  valid  for  m(d), 
any  arbitrary  function  of  Its  argument  6. 

The  Maxwell's  equations  for  this  problem  can  be  expressed  (the  time 
dependence  Is  assumed  to  be  e^)  as 
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7’xE(r)  =  -j<4W<£,(r)  *  H(r)  -  ^  1^ 
yxH(r)  =  ju*Q  fXr)  •  E  (r) 


(2) 


where  the  relative  dielectric  constant  t( r)  and  the  relative  permeability  n(r) 

***  xv* 

are  defined  In  the  following  way. 


^rr 

0 

€(r)  = 

'  Or 

eeo 

0 

0 

0 

‘zz 

,  for  O^r^a 


(3a) 


with  constant  elements 


6  2  (constant)  ,  for  a^  r  ^  b 


rr 


'rO 


'06 


‘Sr 


(3b) 

(3c) 


A*(r)  = 


rr 

* 

<D 

0 

-*er 

0 

,  for  r^a 

(4a) 

0 

0 

Mzz 

with  constant  elements 


(4b) 


I 

, _  THE  UNIVERSITY  OF  MICHIGAN 

4386-1-T 


=  Hi  (constant)  ,  for  a^r<b 


^rr  =  »ee 
^t6  ~  ^ dr 


(4c) 


r  In  eqs.  (2)  represents  a  three  dimensional  position  vector,  and  r  Is  the 
radial  coordinate.  The  results  developed  In  Appendices  A  and  B  will  be  fre¬ 
quently  used  In  the  following. 

A  method  of  solving  any  source  problem  In  terms  of  Green's  function 
will  be  presented  here.  To  construct  a  Green's  function  for  a  problem  with 
some  given  boundary  conditions.  It  Is  sufficient  to  find  corresponding  eigen¬ 
functions  which  form  a  complete  orthogonal  set.  These  eigenfunctions  are  solu¬ 
tions  of  the  source-free  problem  subject  to  the  same  boundary  conditions. 

In  the  present  problem  where  the  waveguide  Is  uniform  (Independent  of  z) 
and  the  medium  Is  also  homogeneous  In  the  axial  direction  z  (l.  e. ,  components 
of  £,  and  £  are  not  functions  of  the  coordinate  z,  with  £  and  /£  having  forms 
shown  In  (3)  and  (4)  respectively),  one  can  assume  that  there  will  be  waves 
propagating  in  the  z-dlrectlon,  having  z-dependence  as  e~^z  ,  where  X  Is  a 
propagation  wave  number  for  a  particular  mode.  This  assumption  leads  Max¬ 
well's  equations  with  appropriate  boundary  conditions,  to  an  eigenvalue  problem, 
wlthX  as  an  eigenvalue  (see  Appendix  B  and  ^9  J  to  [l3]  ). 
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Thus  the  source-free  solutions  (r)  and  (r)  satisfying  the  following 
Maxwell's  equations  (5),  form  a  complete  orthogonal  set  of  eigenfunctions. 


V  x  £jr)  =  -  )W0lAr)  •  Xj^  (r) 

V  x  Xjr)  =  jue0£(r)  •  (l> 


(5  V 


The  orthogonality  relation  can  be  obtained  by  choosing  another  set 
^^/(r)  and  X ^'(r),  which  satisfy  the  same  boundary  conditions  and  the  fol¬ 
lowing  Maxwell's  equations. 


V  x  £^»(r)  =  j<4*0£+  (r)  *  K 
7xX"^(r)  =  -ju€0l+  (r)  '  £  (r) 


(6>+ 


where  *  denotes  complex  conjugate, 

€+  (or  p+)  =  adjoint  of  €  (or  n)  ■  complex  conjugate  of  the  transpose 

/v  ^  /v 

of  e  (or  m). 

a*  X/ 

Now  It  can  be  shown  (see  the  above  mentioned  references)  that  the 
required  orthogonality  relation  Is 


/<£>  *  So  *  *  Ni^  *  jj  • 


(7) 


+In  general  the  single  Index  JL  (ort  )  Is  actually  a  double  Index  In  (or  t'n') 
corresponding  to  radial  and  angular  variations. 
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where  l  and  jt‘  correspond  to  -th  and  X ^  *th  modes  (eigenvalues) 

respectively. 

z  =  unit  vector  In  z -direction 
-o 

=  normalization  constant 

V =  *• for^  =  *** 

=  0,  forX^ 

and  g,  are  transverse  components  of  and  respectively,  s  Is 


:u  — 

the  cross  section  of  the  waveguide. 


Since  In  the  present  problem  the  only  source  Is  a  magnetic  current,  the 
total  fields  can  be  expressed  In  the  following  way,  using  the  appropriate  dyadic 
Green's  functions  (eqs.  (16b)  and  (17a)  In  Appendix  B): 

2ir  b  « 


E(r)  =  - 

and 

H(r)  = 


£;(r)j^/.(rl)  •  gp  Im  (rj)  r'  dz'  dr*  dfl' 


(8) 


2N, 


0  0  -oo 

2  Jr  b  oo 


gjjfejg)  •  gp  1m  (£l)  r'  dz'  dr'  dfl' 

2N/ 


(9) 


0  0  -oo 


r  =  observation  position  vector 

r|_  =  source  position  vector  (l.  e. ,  primed  coordinates  refer  to  source) 
and  the  time  dependence  e^  Is  suppressed  everywhere. 
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Using  eq.  (1),  the  above  two  expressions  can  be  reduced  to  the  following 


forms 


E(r)  =  -cV  L  l(~  \  K  e'  *  O)m(O') 


H(r)  =  -c 


2jt 

5 U>9,-o,m<9'’ 


where  (c,  6' ,  0)  Is  the  ©-component  of  X  j(c,  ©' ,  0).  If  (r),  H 

and  X  £  (r)  can  be  expressed  as  (the  single  Index  JL  Is  replaced  by  the  double 


Index  In): 


llnto  =  Allniln<E> 

.  Ai  [n  Sin  <-> 
K  in <£)  ■  Aln  2  [n  M 


and  moreover  If 


&'ln  <E)  =  £'(n  (r»  z>  e 


l  =  1,  2,  3,  ... 
n  =  0,  il,  tl. 


then  the  above  expressions  (10)  and  (11)  can  be  reduced  further  to  the  following 
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=  "  |  51  f~ln  <£>  «Mn  (c)  ®n 

fi<£>  ■  -|E  (^S~)  «•  ®  *«n  (0)  Kn 


where 


e -JQ0' m(0')  d0'  , 


(16a) 


and  A|^  are  constants  and^  (r)  ,  and  g'^  (r)  are  known  functions. 

If  m(0')  =  m,  a  constant,  then  n  =  0  |  I.  e. ,  —  s  0 1  and  consequently 

1  b6  J 


m  =  2irm 
n 


For  n  =  0,  the  above  summation  Is  no  more  than  a  single  summation  over  l. 
Now  using  (12)  and  the  orthogonality  relation 


Oj  £tln(£)  •  K'ttn  W  *  5o  *  Nto 


It  can  be  shown  easily  that 


A,,  A'' 
lln  lln 


ltln(-  ‘  *  5o  ds  "  2wb 


It  b 

V  j[tr[n(I)*9In(£)',e[a(£)eJln(E)] 
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Now  using  the  above  relation  (18)  In  the  expressions  (14)  and  (15),  It  can  be 
shown  that  the  complete  fields  have  the  following  forms 


E(r) 


H(r) 


c 

2 


c 

2 


z 


_ mn  iln  <£>  gj$ln  (c) _ 

|  ([frto<£,geto(E)  •  rdrd9 

_ Sn4,  <r>  «8[»  (c) _ 

\  fl Tf  (r)g"  (r)  -  f  (r)g"  (r)"|  rdrd© 
}  I L  rln  0ln  0ln  rln  “  J 
''O  J0 


(19) 


(20) 


Although  the  above  two  expressions  (19)  and  (20)  are  valid  for  both  dissipative 
and  non -dissipative  media,  the  relation  between  and  becomes  simple  In 
the  case  of  non -dissipative  medium.  Thus  In  a  non -dissipative  medium  (see 
Appendix  B  and  £  9  J  )  It  can  be  shown  that 

*in  =  tin  =  complex  conjugate  of  gjn 

(21) 

-In  =  -In  =  comP^ex  conjugate  of  fjn 

Therefore,  for  non-dlsslpattve  medium,  the  total  electromagnetic  fields  are 
given  by  the  following  expressions  (using  (21)  ): 
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E(r)  =  -  - 
2 


(22) 


-|E  ITb 


mnKtn(r)  g&ta  (O 


“f 

Jo  Jo 


frto(r)*L(£)  -  «W»  <23) 


It  may  be  noted  here  that  for  non -dissipative  medium  the  propagation 
wave  number  X^n  (eigenvalue)  In  the  z -direction  is  a  real  number. 

An  alternative  set  of  expressions  for  total  electromagnetic  fields  which 
are  particularly  suitable  for  dissipative  media  (although  valid  for  non-dlssipattve 
media  also),  can  be  obtained  from  (19)  and  (20)  using  the  following  transforma¬ 
tions 


(24) 


A  In  =  “  ^ln 

£r0>  "  “ifo  (l*  _^in  “  ^rd’ 

£r«-  *‘re)  5  *ta(£-  A'  -£r#'  *  *ta<£> 

Thus  the  total  fields. which  are  particularly  suitable  for  a  dissipative 
medium,  can  be  expressed  in  the  following  way  (using  (24) ): 
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*0  ■  fZs 


<£>  «9ta (c) 


rdrd9 

Jo  Jo 


v-T  ”,  (r)  8fl«„  (c) 

H(r)  =  — - — - — - 

U  \  f[W^9Cn^'W0*rta®]rdrd8 
Jo  Jo 


(26) 


It  should  be  pointed  out  here  that  for  a  dissipative  medium,  the  propagation 
wave  number  (eigenvalue)  In  the  z-dlrectlon  Is  complex. 


Solutions  of  the  Homogeneous  (source -free)  Maxwell's  Equations 

It  has  been  demonstrated  In  Appendix  A  that  for  source-free  and  homo¬ 
geneous  (or  for  piecewise  constant  n  and  t  )  medium  the  longitudinal  (z-com- 
ponent)  components  of  electric  and  magnetic  fields  obey  the  following  two 
equations 


t  z  «rMr  i  2  CrnT 


a3 


\J  _  ^Z  « 

t  '2  *r**r  &l  z  ~  €r»r  &3  2 


V2V  + 


(27)* 


(28) 


*  To  simplify  notation.  Indices  are  omitted  from  both  field  quantities  and  X  . 
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where  =  k26r  (n ^  -^,2)  -Mr^2 

ai  =  k%<fcr  '£'2)  "€/2 

a3  *  i«r€'  +*i'£r 


Coo  a  C, 


^rr  ■  €ee 


(r$  3  *8r  =  fc' 

“  *2 


“rr  *  “««  3  “r 


"rO  3  “Or  = 


fi  -  jU 
ZZ  2 

,2  2 
k  ~u/je 
o  o 


The  transverse  fields  (l.e. ,  r  and  6  components)  can  now  be  expressed 
(see  Appendix  A)  In  terms  of  <£z  and  Xz  In  the  following  manner 

4,  3  £[*MA -  U9) 


ir*  ■  )utoai'7t^-l  (30) 
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=  k2C'  -  a*’2)  +X.2m* 

a,  *  k2^'(€2  -€  ’2)  +X2*' 

Z  r 

a.  =  +#i'€')  -X2 

4  r  r 

Pi  =  k4a§  -  a| 

It  Is  observed  here  that  If  =  0,  the  two  equations  (27)  and  (28)  become 

uncoupled.  *  Although  this  Is  a  necessary  condition  that  the  conventional  E-type 

and  H-type  modes  separate,  It  Is  not  sufficient.  The  sufficient  condition  depends 

on  the  boundary  conditions.  For  example.  If  the  medium  completely  fills  a 

perfectly  conducting  waveguide,  and  If  =  0,  E-type  and  H-type  modes  can 

exist  separately.  But  on  the  other  hand.  If  there  are  two  coaxial  media,  (the 

outer  one  may  or  may  not  be  bounded  by  a  perfect  conductor)  E-type  and  H-type 

modes  can  exist  separately  If  and  only  If  ~rr~  =  0  (with  a3  =  0).  But  If 

ou  ° 

4  0,  and  even  If  a^  =  0,  In  the  above  two  coaxial  media -system  E-type 

and  H-type  modes  cannot  exist  separately.  A  similar  discussion  for  Isotropic 

media  where  =  0,  can  be  found  In  fl6j,  sec.  11.6. 

For  the  solution  of  (27)  and  (28)  £ z  and  Xz  can  be  eliminated  yielding 

a  single  4th  degree  equation  In  each  of  t  and  V  which  satisfy  both  of  die 

z  z 

second  degree  equations  (27)  and  (28).  If  a  choice  Is  made  such  that0=£  , 

z  z 

*  TE  and  TM  modes  also  decouple,  l.e. ,  they  can  exist  separately  if  there 
Is  a  constant  line  source  In  the  z-dlrectlon.  In  this  case  it  =  0. 
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then  the  equations  (27)  and  (28)  can  be  reduced  to  the  following  equation  (for 


detail  see  Appendix  A): 


V2  0  +  yfz  0 


where 


11  -“VVJ  *  r  !vf  L 1 =  ~ 5 ~  J 


Solving  equation  (32)  for  a,  one  obtains 


€ z al  ~**zal  +  [(al ValMz)2  +  ^ zj 


2u  (  €  a,/ 

O  Z  o 


Therefore,  the  roots  of  ip  can  also  be  expressed  In  the  following  way 


yj'  =  v  +  1/v2  -  u 

*1.3 


where 


»1  «.♦«!<*»  = 

2  £  «r  "  2 
r  r 


(35a) 


«z*z 

crMr 


-k4,l] 


€ZMZ  »2  ,2 

Mr‘r  =  7  a 
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where  J  and  N  are  Bessel's  functions  of  the  1  st  kind  and  2nd  kind  of  order  n 
n  n 

respectively. 


Since  0!  =  6z  +  jajVz  =  Jn  tyi  r)eJn®,  and 

0 s  =  t +  j  v„ =  ^  r) e  * 


It  Is  easy  to  verify  that 


£z  =  - r)  -  aj  A!  Jn(?fi  r)J,  for  O^r^a, 


(41) 


and 


.  .  jn0  r  1 

Hz  =  [A*  Jn  (7a  r)  -  A!  Jn  (7f  r)J  ,  for  0^  r  ^.a. 


(42) 


Since  the  boundary  condition  requires  that 


£  (r) 

z 


-  °  -  ir  X (r)/ 

or  25 


r=b 


r=b 


the  equations  (39)  and  (40)  can  be  rewritten  In  the  following  way 


*  B,  yn(r)el“9. 

for  a$  r^b 

(43) 

Hz  =  B,  G„  (r)  eJ"S  , 

for  a^  r  £  b 

(44) 

where 
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/n(r)  =  Jn  W  Nn  <*r)  -  Jn  (t}r)  NQ  (r,b) 

(45a) 

Gn<r>  "  Jn  (*r)  (rlb)  -  Jn  (,Tb)  Nn  <7Ir) 

(45b) 

B,  =  -  B'" 

Nq  (»;b) 

(45c) 

BL" 

B»  =  — 3 — 

N'  (>fl>) 
n 

(45d) 

/  dNn  (Jfr) 

-  <$r 

r=b 

(45e) 

j.  (,bi .  d,"(’»r) 

n  1  d  (7)r) 

r=b 

(45f) 

Now  using  equations  (41)  to  (44),  In  the  relations  (29)  and  (30),  the  transverse 
components  of  the  source  free  solutions  of  Maxwell's  equation  (5)  can  be  expressed 
as 

..J“®  r  f  .  .T_  {*:  r)~\ 


'r  Pi  («!-«,)  €. 


H,  ,  J„  (rjl  r)~l 

€zTi  R  W  r>  -  VV?2  — ; —  I 

f  ,  Jn  1*1*  r)  I 

►  A,  |  Vz  T  Ja  (Vi  T>  +  “VrVl  - r - j j  ’ 


for  04  r<  a 


-  THE  UNIVERSITY  OF  MICHIGAN 

4386-1-T 


r°^3  B,  e^n0  Gn(r)  - Bi  e^n0  Cn(r),  fora^r^b 


Jn9  f  f  ,  j  (ift  r)- 

^^|-^r€r^2  Jn^l  r)+**zR  ' 


+  aJ  UM0*tMt7l'j[  i\ '2Jn%  r)  +  ntzT  - —2—  |  for  04  a 


_(n'  r)T 

— —2 —  for  04  r  4 
r 

-J  a 


V  ®2  ®  sn  <r)+-afL  Bi  e 

7 


jn0  */n^ 


,  for  a^  r£  b 


(«!  -<*3) 


^  R' 


J'tyl  r)+ 


n^2  Mtr  *4  Jntyi  r> 
*  €z  a3  r 


f  ,  1  n7l  Str*4  Jn^a 

A3  <  ^3  T'  JjJ  (73  r)+  ~~z -  ’  for  0^  r^, 

L  *  €z  ^2  a3 r  J  J 


nu€0€j 


Bj  e^Vn  (r)  + 


j^ein6 


for  04  r  <.  a  (48a) 


SQ  (r)  ,  for  a^  r^  b  (48b) 


i£-U 


3  3  f  *51 

€r^r^l72*1*  »  nR 


r>} 


for  04  r4  a 


(49a) 


!2¥ii.«M0LW*4t^^,«Urt. 


for  a.<  r^b 
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where 


“r£r 


M  = 


Czai 

Mr€r 


*2 


R  =  a4<*3  -umq  aj  = 


T  =  ®3  -/a4«i 


R’  =  ue.aa^-^a. 


T'  =  uCo  = 


a4VrM  "^c2a3«« 
wcocza3 

= 

w«o«za3 

M  aaMrCr  -^a3a4ez 
*  a3*z 

®  a3a4^z 

3  z 


(50a) 

(50b) 

(50c) 

(50d) 

(50e) 

(50f) 


sn(r)=  J;(7b)^(^r)-J»1(^r)N'l(7b)  =  -  j  ~  Gn(r) 

Cn(r)=  Jn(77b)Nn(7r)"Jn(7r)Nn(7b)  3  y  ^  / 1(r) 

Dispersion  Relation 

Since  In  the  present  problem  the  eigenvalues^  are  discrete,  the  boundary 
conditions  satisfied  by  the  total  fields  (l.  e. ,  the  fields  due  to  the  presence  of 
the  source)  are  the  same  as  those  satisfied  by  any  Individual  fields  (l.  e. ,  source- 
free  field).  Therefore,  the  following  boundary  conditions  can  be  Imposed  upon 


(50g) 

(50h) 
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the  source-free  fields,  6Z,  M z ,  and 


6z(.-)  -  Szuh 

at  r  =  a 

(51a) 

60(a -)  =  6e(a+) 

at  r  =  a 

(51b) 

£z(b)  =  0 

at  r  =  b 

(51c) 

(b)  =  0 

at  r  =  b 

(51  d) 

U'i  (a*)  =^z(a+) 

at  r  =  a 

(51e) 

^(a-)  =^0(a^) 

at  r  =  a 

(51f) 

The  constructions  of  and  €-q  In  a^  r^.b  are  made  In  such  a  way  that 
the  boundary  conditions  (51c)  and  (51d)  are  now  automatically  satisfied,  since 
.^^(b^OS^b).  If  the  remaining  boundary  conditions  In  (51)  are  Imposed  upon 
the  source-free  fields  expressed  In  the  equations  (41)  to  (44)  and  (47)  and  (49), 
the  following  relations  among  the  arbitrary  coefficients  At,  Aj,  ^  and  are 
obtained: 


ai  AaJntyi  al-ofcAjJnty 


i  a)  j  =  Bj^j(a), 

1  f  r  ,2  JnWi  a>  1 

Z^Tz\  Ai  a)+nczR  — - j 

[  i  a)  1  1 

WMo^r^r^i^a  Jhtf a  *)+ncz  T  ^  J  I 


“J^o^a  njf  Vn(a) 

— a,sn(.)+-^.B1 


a 


(52a) 

(52b) 
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«l  -<*2 


A3  Jn0?3  a)  "  Ai  Jn^i 


=  BjGn(a), 


-J 


-Jw£0^a  njf 

=  - W — B1Cn(a)  +  -p-  B,Gn(a)  .  (52d) 

•  y)  a 

Non-trlvlal  solutions  for  the  constants  Ax,  Aj,  B1,  and  B|  exist.  If  and 
only  If  the  determinant  of  the  coefficients  of  these  constants  appearing  In  equations 
(52)  vanishes.  The  vanishing  condition  of  the  determinant  gives  the  characteristic 
equation  (or  dispersion  relation).  Instead  of  calculating  the  determinant  of  the 
coefficients  and  then  equating  It  to  zero,  one  can  also  eliminate  Bj  and  from 
(52)  and  from  the  remaining  equations  It  Is  easy  to  obtain  two  Independent  values 
of  the  ratio  Aj/Aj.  Now  equating  these  two  values  of  Aj/Aj,  one  obtains  the 
desired  dispersion  relation  expressed  In  the  following: 
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*  When  =  1  =  fc3,  the  dispersion  relation  (54)  agrees  with  that  of  Van  Trier 
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It  should  be  noted  here  that  to  obtain  the  expressions  (53)  and  (54),  some 
useful  relations  (tabulated  In  Appendix  A)  have  been  used.  These  relations  will 
be  used  frequently  In  subsequent  derivations. 

When  the  ring-source  m(0)  Is  constant,  n=0  (l.e. ,  -^-EO).  In  this  case 
It  can  be  shown  that  the  dispersion  relation  simplifies  to 


S{^z  Gb(a)Jifai  aHi^2So(a)Jo(Tha)}  _  i(*)ia)+V7l  CQ(a)J0(>ha) 


M(yz  G0(a)J1(Tjj  a)-^3S0(a)Jo^a)}  JJ€z/0(a)  Jx  (ijjaHCjjjj  C0(a)J0(^3a) 


rr  (55) 


Alternatively,  equation  (55)  can  also  be  written  In  the  following  manner: 


-  Vt1-  G0(a)/(a)J1(t>ia)J1(^a) C0(a)S0(a)J0(ijia)J0(^a) 
IL  u  w  £  u 


aWa2'  H) 


JQ(^la)Ji07aa) 


J0^2a)Jl^la)  r  /I 

+  - [€3^ZSCo(a)Go(a)+^«ZM^(a)S0(a)J  =  0 


MzMC0(a)G0(a)+Ma<z 


!^(s)SoWj 


A  number  of  dispersion  relations  for  various  special  cases  has  been  developed 
In  Appendix  C. 

The  solution  of  the  dispersion  relation  (54)  together  with  the  relations  (34) 

and  (37b)  gives  an  Infinite  number  of  discrete  values  of  rf[ ,  rj'2  ,  *}  (and  hence  X 

also).  The  radial  wave  numbers  rj  J  and  *)\  can  also  be  called  eigenvalues  In  the 

domain  0£r^a,  and  similarly  rj  is  the  eigenvalue  in  the  domain  a  <r^_b,  of  the 
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where  either  the  value  (57)  or  (58)  may  be  assigned  to?la  since  they  are  equivalent  through  the  dispersion  relation 
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Now  using  (59)  and  the  relations  (52a),  (52c)  of  this  chapter  and  (38(4)), 
(38(10)),  and  (38(11) )  of  Appendix  A,  one  can  express  Bx  and  B}  In  the  following 
manner: 


Bi  =  At 

B,  =  Ax?3 


where 


S  ti  a)  -  M  Jn(iji  a) 
(7a  -7i3)  </(a) 


a)  -  JQ 


(60) 

(61) 


(62) 


(63) 


Expressions  for  Source-Free  Fields  In  Terms  of  Only  One  Unknown  Constant  Ai 
Since  all  the  unknown  coefficients  A3,  Bj  and  Bj  are  now  expressible  In 
terms  of  the  only  one  unknown  Ax,  the  source-free  fields  can  be  written  In  the 
following  way  (e  ~^z_wt^ls  assumed  to  multiply  all  the  expressions  for  the  fields): 

AieJna  r  ,  1 

t  =  — -  S?x  Jnrt3r)  -  MJn(f,r)  =  A^,  forO^r^a  (64a) 

<7a-7ia>  L 


Aj  fz,  for  a^  r^  b, 


(64b) 


32 


THE  UNIVERSITY  OF  MICHIGAN 

4386-1 -T 


tT  = 

r  a  a 


TT~S - 7~  {  *zl  1 B  Jn<7‘ r)  ‘  ”ul*'f-rUI’) a  ——} 

i7a  ^a  ">7i  )  L*-  r  J 


*r^r7i7a  tya  “7i 


(65a) 


+  ?ijVa  T «zJ'n(7i r)  +tK4i0€I^r7i2  -Jp(^  )  jj  =  Aifpt  ...  for  0$  r£a 
f  ni4ijx2  1 

6r  =  Aie  U  jj-y -  "  ~  %2  Cn^rM  =  Aifr»  •  •  •  for  a£  r^b  (65 


*6=  ~ 


■  z^zu<p^a3e^P9Al 


2  a  >2  i3 .  |3  i2.  n 
cr  **r7i  7a  tya  ~7j  )  It 


+  n€zR 


Jntyl r) 


} 


+  ^i|^<4^Vi*7i372  Jnfyar)  +  n€zT  D  ^  =  Axf0  ,  for  O^r^a, 

c  jn©[  ~j^c^2  «  n/?,  /  1 

*0  =  Aie  ^ -  *3Sn^  +  — a - ^n^rM  =  Axf0,  for  a*r<.b,  (66b) 


Jw€o/a3€zAieJne 


^rer<7a2  -?*> 


-  Jn(7lr) 


Ki  Jn^Ta  r)  “  Jn(7ir)  -  Axg  ,  for  0£r£a,  (67a) 


*4  =  ai  VJD0Gn(r)  =  Algz*  f°ra^r6b. 


(67b) 
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"r **>*)•  ~  J{*‘E*V‘»»3Jn<'>ir)  +  “s~'| 

«r  »rli  7i  Wi  -7.  i  J 


(68a) 


-  fijijTVtej  +  “7iJ3‘J“r  ~ST~}  ■  Ai*r  •  for  0  -r£a- 


*r  =  Aleine  ^  ?, /(r>  +  Sn(r)l  =  A,gr, 

7  r  7  J 


jn0 


fora^r^b,  (68b) 


_  i^o^z  Ai e  |~p  *  ,  Jn(7lr) 

*  "  ,»  *v  H  MJn<7ir)  +  QR^£z^za3 - 1 - 

£rMr7i72<72  "7i  '  lL  r 

^  f  9  9  (9  l  f  A  ^  }1  I 

-  « i  |v*r9t  72  s  Jn(72r)  +  nT'J(eznza3  — -  JJ  =  Aj  gg,  for  0  4  r 


} 


(69a) 

4a 


*&■  A 


1n6  I  "jufeo£2  %  ntf  ¥ 

i eJ  - - - Cn(r)  +  — j — —  G,/r)  =  Aigg,  for  a^.  r  4b,  (69b) 

L  v  T  r  J 


Determination  of  the  Constant  At 

Determination  of  the  unknown  constant  Ax  depends  on  the  orthogonality 

condition  satisfied  by  the  source -free  fields.  Since  orthogonality  conditions  are 

different*  for  a  dissipative  and  a  non -dissipative  medium,  the  constant  Ax  will  also 

*  Although  the  forms  of  the  orthogonality  relation  given  In  (22)  of  Appendix  B  are 
valid  for  both  dissipative  and  non-dlsstpatlve  media,  they  are  particularly 
suitable  for  dissipative  media. 
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be  different  for  the  above-mentioned  two  media.  Here  Ax  will  be  calculated  for 
non-dlsslpatlve  medium  only.  For  dissipative  media  the  corresponding  A!  can  be 
calculated  using  equation  (18)  together  with  the  transformations  given  In  (24).  The 
primary  aim  here  Is  to  calculate  the  total  fields  with  amplitudes  due  to  a  given 


source.  Therefore,  to  calculate  total  fields  It  Is  only  necessary  to  find  the  ratio 


Nin 


for  any  mode  (eigenvalue)  In.  This  ratio  can  be  determined  from  the 


knowledge  of  source-free  fields,  as  suggested  In  equation  (18)  together  with  (21). 

aJ1 


In  other  words,  the  ratio 

i  12 


N, 


In 


Is  given  by 


1  In 


N 


In 


r2n  A) 

\  d0\  [frln(£>  &Wl)J  r  ^ 


(70) 


The  above  Integral  In  (70)  can  be  expressed  in  the  following  compact  form 


32 


Nln  F 

,  ta  =  \  < 

l8(  r,rto<£)«%Ui®-,eto<£>«rlii<^] 

o 

“¥ 
T— 1 
< 

)o 

:=i 


(71) 


where 


Ftnl  =  LlnV*rl  ’)1lnf€z[cr^r7alnMin^n  "  ^0^2inRin  (72-1) 

Fin2  =  ^Ln^r^rliLn^Jln  tilntz[^o7aln’^ln^lif<za3  “  ^In^i^rfttn^lnJ  (72-2) 
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Fta3  '  LtaW.?ito7ito».lnLTto^^to«to+*oA“lrfV2*3]  <72-3) 


[n4 


Ln5 


ln6 


ln7 


ln8 


ln9 


VA*^ «  Jl^ln  fE^ol  + Cr“rTln7itaSln]  (72-4) 

l?3U.flNn(Jlnb»|1]  <72‘5> 

TT^  [•.-Hj  W’ VV>4V»  |*>lnf 

I  “I  (72-6) 


U^. 

pin 


fro'ft.f  VV>W,4V»  |?3ln|J£n(7tob)V7tob»3 

(72-7) 


j^z  foetal  l,n<V>,l  +'‘^|<3ta|K«lnbfl 


(72-8) 


*  nLin7il„  -  “<‘o‘‘rtrl»’iJ’  “to]  <72-9» 


F.  =  nL.  y\ ! .  |  fu»M  e  V  V) 12  i[2.  S{  -€?R.  T/nXL/u  a„"|  (72-10) 

mlO  in'  in  inL  o  r  r  '2in  l1  in  in  z  in  in  in  z  3J 


.3  3  »2  * 

\» 


Flnll  =  ”Lln1u„&z«toRin^z*3-“*‘o«?“r|7iJi«r  fiin  ]  m-ll) 


Ftol2  =  ”  LinV‘rtz[Rln  C  “inVr  ’  "M»  to  "Xf»*s]  <72‘12> 

*«*r  3  3  *  *2  i2  *  2  *  t  *  “1 

Finl3  =  ^inJ^n^inL^oS^rViin^Jin^  “<zRinTii^in#i2a3j  (72-13) 


t2  t 
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*  %  r  *2  *»  -  *t2  *  - 

"  VrW'taK’J’taJ'ta  <72-14 

nLin1lln5‘lnC“,‘o<r'‘??il,lniataMin  +  TinRto*^^  J  (72-15) 

n^ln^aiiJ^1inj*tcz7,in7'lri^in*i2a3  +  tJ*lo*r^rH1  Inf  ®lnj  (72-16) 

“LJ‘lnf>ln[«2  Tin&(Xi'‘z,3  +  “Vr^rl'4,  J  <72-17> 

2  P"  *  *2  ♦  12  jflc  f  mm 

“  V‘  inCzt‘rSLTln,lau1Miner  "r  +  “w<?i  ^ln^lnVaJ  <72-W 

J  [‘iTlnTL,^nMI«3  +  u“oelMr  jl’hf  ^ta]  (72-19) 

Li„  | 5  J’  £r,'r£l["Wl„K,'.*3  +  TU>?>toWr]  <72-2»> 

^  Ck,^*,53ta?21„Jn(V)^<,llnb)-45anf31n3r(,Ilnb>Nn'V)! 

(72-21) 

3U>V \b)V  *3  J2i»Mn<  VA*  7mb)] 

(72-22) 
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?in24  ■  ^[k^‘»53ln\in 


(72-24) 


in25 


In  *in  /-„  , 


(72-25) 


ln26  =  -^T-  Ck2^£>53j2inV7u.b)V\.b,-^2  Ja^toWN^b)] 


(72-26) 


■ta27  *  T%^k2b»<>53J21nJn(11„bK<V><52j31nJn<V>VV>] 
win  , 


(72-27) 


jn  _  2  *  *' 


*  * 


in28 


fHfin52ta531nVC^(V,-k^»531I.52U,VV,Nn<V»] 

in  in  /7o_<« 


in29 


in30 


in31 


in32 


(72-28) 

(72-29) 

(72-30) 


M>.V>f + ^|f3to|V>tob>f  ] 

2  ^ 

wf*  ^3j\wn 

2 

nr  [‘o‘»k|V Vh iW W1 wO 


'in 


(72-31) 

>] 

(72-32) 
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(74-27) 


(74-28) 


(74-29) 


(74-30) 


(74-31) 


(74-32) 


It  may  be  noted  here  that  for  n  =  0  (i.  e. ,  when  —  =  0),  F.^  =  0,  for  JL  )  9, 
and  all  the  Integrals  for  X  £  8  can  be  evaluated  in  closed  form.  But  when 
n  ^  0,  the  above  integrals  can  be  expressed  partly  in  closed  form,  partly  in 
series,  or  they  may  be  calculated  numerically  also. 
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Expressions  for  Total  Fields  (for  loes-less  media) 

Now  substituting  the  expressions  (64)  to  (69)  and  (71),  in  the  equations 
(22)  and  (23),  the  complete  fields  due  to  the  magnetic  current  ring  source  can  be 
expressed  In  the  following  way  (suppressing  the  time  dependence  factor  e^: 


_c_  V"  — 1 
'  4*  4i 


+  )■>» 

m  e 
n  _ 


1,  n  Vi  FinAn l 


[s,  Si.  J  tya.  r)  -  M  J  (Vl4  r)  1- 

in  *in  n  ,ain  In  n  ,lln  J 


(75a) 


( 


i“£oc^iflc*m(c)  ,  **J3Ji*M 

+  - IT* - 


Tin 


,  forO<r<a 


E 


m  e 
n 


+  j/iQZ  +  jnd 


FlniIinZ' 


?V^inr)- 


(75b) 


in 


- j? -  .  for  a*  r4b 

cl?ln  J 
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*1in^ain^aln"^lin) 


JL 


Vnv> 


U  *U 

t.  - -  Ufi  n 

in  r  o 


r£r^ta'!’inJn(,fllnr)) 


-  ,7  r>f  n  e  T  i&£}l 

r  r  'iin  2i“  n  2in  z  in  r  J  I  * 


for  04  r<  a 


+»V*»  r 

srTTL 

'  lnttrul 


1UCoC^,nClnfe) 

JL  T 


>W1 

*3 

C  It  J 


liilA1.  JU“oMi5310Sln<r) 


fl 


for  &At< b 
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HIGA 


*  * 
^3to  ■ 


*V,  tS*  V^AnS^int'U 

l*i 


°*J3iA><c>l  fr> . . . 

“If - J  QVr1«I«%A.VV>+l*b 


i  J  (ft.  r)  -| 

iP  *  n  4n  l 

in  in  ^ za3  J 

T»  -» 


-  ^VXVt.'.^A—  Vn<r>l  I 


for  O^rAa 
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- - r  m  e 

H  *  -  — \ 

9  4*  / 

XT 


z+jn0  *  *  . 

iU(oC‘%inCinM  , 

- ; - + - p - 

•  v  <  ■ 


x  +  Ayv?  1 

L  i..  v  J 


(80b) 


Wherever  the  sign  ±  or  +  appears,  the  upper  sign  corresponds  to  the  propagation 
in  the  positive  z-direction  and  the  lower  sign  for  the  negative  z-direction. 

Expression  for  Average  Power- Flow  Due  to  a  Magnetic  Current  Ring  Source 
The  average  power  flow  is  defined  as 


P  =  ±  Re 
av  2 


51  i>i 


2  dS 
o 


x  z  dS 
~o 


(81) 


where  Re  means  real  part  of 

E_  «  total  electric  field  due  to  the  source  at  any  point 
H.  =  total  magnetic  field  due  to  the  source  at  any  point 

Now  using  equations  (18),  (21),  (22)  and  (23)  in  (81),  the  expression  for 

P  can  be  written  in  the  following  manner 
av 
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In  particular  when  n  =  0,  i.  e.,  when  the  ring  source  is  of  constant 
amplitude,  the  expression  for  the  power  flow  reduces  to  the  following  form: 


“V*  Cio(c) 

23  FiotlioJi 
1*1 


It  may  be  noted  that  for  non-dissipative  media  the  quantity  inside  the 
square  bracket  is  real.  Any  individual  term  in  the  series  in  (82)  or  (83) 
represents  average  power  flow  corresponding  to  that  particular  mode  in  (or  i 
when  n  =  0). 
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n 

WAVE  PROPAGATION  IN  AN  ANISOTROPIC  PLASMA: 
SLOW  SURFACE  WAVES 


Introduction 

In  this  chapter  the  general  results  of  the  previous  chapter  will  be  applied 
to  the  study  of  propagation  of  electromagnetic  waves  in  an  infinitely  long  aniso¬ 
tropic  plasma  column  enclosed  by  a  dielectric  cylinder,  which  is  also  enclosed  by 
a  perfectly  conducting  metallic  cylindrical  waveguide;  i.  e.  the  geometry  and  the 
source  of  excitation  are  the  same  as  those  of  the  general  probelm  except  that  in 
the  present  situation  the  anisotropic  medium  is  represented  by  a  plasma  column 
with  a  uniform  static  magnetic  field  in  the  axial  direction  z.  The  relative  permea¬ 
bility  (I,  of  the  dielectric  medium  which  encloses  the  plasma  column  is  assumed 
to  be  unity. 

The  plasma  is  considered  to  be  fully  ionized  (i.  e.  macroscopically  neutral) 
and  there  is  no  drift  velocity  (d.  •*  )  of  electrons  or  of  ions,  i.  e.  the  plasma  is  also 
stationary.  If  one  also  assumes  ...at  the  illuminating  electromagnetic  waves  are 
weak,  then  it  is  possible  to  describe  a  plasma  as  a  dielectric  medium.  In  this 
analysis  it  will  be  assumed  that  the  plasma  is  homogeneous,  i.  e.  its  density  (and 
hence  dielectric  constant)  is  not  a  function  of  space. 
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In  the  presence  of  a  static  magnetic  field,  the  dielectric  constant  of  the 
plasma  becomes  a  tensor,  which  means  it  is  an  anisotropic  medium,  if  the  static 


magnetic  field  is  applied  in  the  axial  direction,  it  can  be  shown  2  3  1 5  that  the 


plasma  has  the  following  dielectric  tensor 


e  je 

rr  J  i 


^  -J'er  £« 


where 


err  =  *dd  =  cr  "  1  + 


2/,  J*\ 
J-Jv-  j*-)2 


uc  UP 

€r*  =  eflr  =  e' =  —  •  T~T7~  Tv 

u  U)  -  U  (1 - ) 

C  U) 


ezz  "  €z  ‘  1  ■  2  n  iv 

u  (1 - ) 

u 


v  -  collision  frequency  (radian) 

% 

- - =  cyclotron  frequency  (radian) 

c  me 

qe  =  charge  of  an  electron 


mg  =  mass  of  an  electron. 
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B0  =  d.  c.  magnetic  induction 


(3e) 


/  q"  n  \  1/2 

u)  =  I—® — “-]  =  electron-plasma  frequency  (radian)  (3f) 

P  \eo  “e/ 


Ng  =  electron  density 


(3g) 

eQ  =  free -space  dielectric  constant  (3h) 

In  the  above  analysis  the  motion  of  an  ion  due  to  a  disturbance  is  neglected 
in  comparison  to  that  of  an  electron.  The  relative  permeability  of  the  plasma  is 
assumed  to  be  unity. 

It  can  be  shown  from  the  relations  in  (2)  that  the  components  of  e  satisfy 
the  following  relation 

e'2  =  (1  -  er)  (ez  -  er)  (4) 


An  interesting  conclusion  can  be  made  from  the  relation  (4),  namely  e'  =  0,  for 

either  cz  =  er  or  cr  =  The  physical  interpretation  of  these  results  can  be  given 

in  the  following  way.  For  isotropic  plasma  (i.  e.  when  Brt  =  0),  e  =  e  ,  and  e'  =  0. 

o  z  r 

On  the  other  hand  when  B^  -f  oo  (i.  e.  u  -*  oo),  c  -*1  and  e'  0.  The  above 

o  c  r 

statements  can  also  be  verified  directly  from  (2).  ft  may  be  noted  here  that  the 
collision  term  v  in  the  expression  for  e  in  (2),  represents  loss  in  the  plasma. 
Although  the  various  dispersion  relations  developed  in  Chapter  I  and  in  Appendix 
C  are  valid  for  a  lossy  plasma,  the  complete  field  expressions  obtained  in  the 
previous  chapter  are  not.  On  the  other  hand  if  it  is  desirable  to  find  the  expres¬ 
sions  for  a  dissipative  medium  (i.  e.  plasma),  one  must  use  the  appropriate 
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orthogonality  condition  and  the  resulting  Green's  functions  which  have  been  dis¬ 
cussed  in  Chapter  I  as  well  as  in  Appendix  B.  But  here  only  the  loss-free  ( v  =  0) 
plasma  will  be  considered. 

Although  the  results  obtained  in  the  previous  chapter  are  valid  for  all 
possible  modes  of  propagation  in  the  structure,  in  the  present  chapter  attention 
will  be  directed  to  the  analysis  connected  with  slow  surface  waves.  The  slow 
surface  waves  are  those  waves  which  decay  radially  in  the  dielectric  region  and 
propagate  along  the  interface  of  the  plasma  column  and  the  dielectric,  in  other 


words  for  such  slow  waves  one  finds 


k 


>  1,  where  Cj  is  the  relative  di¬ 


electric  constant  of  the  medium  surrounding  the  plasma  column.  This  medium 
may  represent  a  glass  tube.  Various  passbands  for  slow-wave  propagation  will 
be  obtained  in  the  following  investigation.  Some  of  the  passbands  depend  on  the 


range  of  the  values  of  the  ratio 


k  /cj 


>  1,  and  some  passbands  do  not  depend  on 


the  particular  values  of  this  ratio,  provided 


k 


>  1,  the  condition  for  the 


existence  of  slow  surface  waves. 

Finally  numerical  computation  will  be  made  for  a  special  case. 


Conditions  for  Slow-Wave  Propagation  and  Determination  of  Various  Passbands 

In  the  following  analysis  only  the  expressions  for  the  radial  wave  numbers, 

rj'  and  rf ,  will  be  considered.  These  quantities  do  not  depend  explicitly  on  any 
1  2 

particular  boundary,  except  that  the  geometry  is  cylindrical  and  uniform  in  the 
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z-direction,  thus  the  results  will  be  true  for  any  such  structure,  closed  or  open, 
which  can  support  slow  waves.  For  a  particular  structure,  one  must  consider 
the  solutions  of  the  respective  dispersion  relation  together  with  the  following 
analysis.  Since  the  solution  of  any  dispersion  relation,  determines  only  a  particu¬ 
lar  set  of  values  of  n*  and  r^,  the  limitation  imposed  on  the  values  of  q|  and  qj, 
which  are  obtained  from  the  expressions  for  q| 2  and  qj2  (valid  for  unbounded 
medium  also)  is  furthermore  narrowed.  Therefore  the  requirements  which  are 
obtained  from  a  study  of  the  expressions  for  qj  2  and  q'2  alone,  are  nothing  but 
necessary  conditions  of  wave  propagation.  The  sufficient  conditions  for  propaga¬ 
tion  of  waves  are  provided  only  by  the  simultaneous  solutions  of  the  respective 
dispersion  relation  and  the  expressions  for  q|2  and  qj2 . 

The  expressions  for  radial  propagation  wave  numbers  q|  and  qj  can  be 
written  here  in  the  following  way 


q'2  =  v 


t 


v  -  u 


(5a) 


k2(l-y){(/32-l)(l-x)+y}  k2xy(/32-l) 


x  +  y  -  1 


2(x  +  y  -  1) 


(5b) 


+  k2y 


(j32-l)2x2  +  4/32x(l-y) 


2(x  +  y  -  1 ) 


where 
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v  = 


2(x+y  - 1) 


|^02  - 1)  -{xy  -2(x+y  -  1)}-  +  2y(l  -y)J 


u  =  - 


k4<1-y>  r„2  „2 


x+  y  -  1 
2 


I 


(/3  -1)  (1 -x)  +  2y(/3  -1)  +  y 


wc 


2 

"p 

y~U 


e  =  1  -  y 


x  +  y  -  1 


c  - - 

r  x  -  1 


,,  ,  j£L 


x  -  1 


/3  =  &- 

k 


In  the  above  definitions  of  ez,  er,  and  e',  the  collision  frequency  v 
neglected. 

It  can  also  be  shown  that 


v2  -  u  = 


.2  1 

2 

k  y 

• 

032  -  l)2  x2  +  402  x(l  -y) 

_2(x+  y  -  1)  . 

In  the  dielectric  region 


(6a) 

(6b) 

(7a) 

(7b) 

(7c) 

(7d) 

(7e) 

(7f) 

is 

(8) 
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n2  =  k2  €,  -  *2 


since 

p2  =  1. 

2 

Therefore  for  slow  surface  waves,  r\  <0, 

2 

& 

I 

II 

1 

k2  €3 

i.  e. 

at2 

n 

-  >  1 

k  ea 

or 

11 

02>  e2 

where 

cs 

II 

«_L 

Or* 

,  S  >  0. 

-  1  >  0 


(10) 


ft  will  also  be  assumed  that  ^  >  1. 

Since  it  is  assumed  that  both  plasma  and  dielectric  are  non -dissipative,  the 

propagation  wave  number  is  always  real.  Although  the  expressions  for  r/2  and 

rf2  show  that  these  radial  wave  numbers  may  be  complex,  on  physical  grounds 
2 

only  the  real  values  of  17J2  and  r£2  will  be  allowed.  For  example*  if  complex 
values  of  tiJ2  and  17J2  (which  are  complex  conjugates  of  one  another)  are  allowed, 
this  means  that  there  exist  growing  waves  showing  instability  of  the  plasma  in  the 


*  The  primary  reason  for  allowing  only  the  real  values  of  n*2  and  Tj2  in  a  non- 
dissipative  medium,  is  that  the  power  flow,  the  characteristic  impedance,  etc. 
must  be  real  for  such  a  medium. 
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radial  direction.  Therefore  this  apparent  inconsistent  situation  will  be  avoided 
by  allowing  only  the  real  values  of  nj2  and  r£a  in  the  following  analysis. 

The  following  three  cases  specify  the  conditions  for  which  rjj3  and  i?’a  are 


real. 
Case  1 


Case  n 


Case  111 


and 


and 


u  <  0 


u  >  0 
v  >  0, 
v^  -u  >  0 


u  >  0, 
v  <  0, 
v2  -  u  >  0 


in  this  case 

rj[a  >  0  and 

r}'3  <  0. 

2 

in  this  case 
>7[3  >  0,  and 

Tjf  >  0. 


in  this  case 

rj[2  <  0,  and 

rfi  <  0. 

2 


(11a) 


(lib) 


(11c) 


The  regions  of  u  and  v  in  the  above  three  cases  are  shown  in  Figure  2. 
r?J2  and  r£3  become  equal  on  the  parabola  v^  =  u  and  both  of  them  assume  complex 
values  inside  parabola,  which  contains  the  positive  axis  of  u.  Therefore  this  is 
the  forbidden  zone  for  u  and  v.  When  y  <  1  (i.  e.  when  <  u),  it  is  easy  to  show 
from  (5)  that  rj[3  and  r£3  are  both  real.  So  even  if  one  considers  instability  in 
plasma,  this  does  not  occur  for  up  <  u,  provided  the  medium  is  loss-free. 
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In  the  following  detailed  analysis  the  possibility  of  the  above  three 

individual  cases,  i.  e.  the  conditions  under  which  the  parameters  x,  y,  0,  etc. 

have  to  be  chosen,  will  be  shown.  These  conditions  on  the  parameters,  which 

are  necessary,  will  give  the  maximum  passbands  of  slow  wave  propagation. 

Before  starting  the  actual  analysis,  it  will  be  convenient  to  introduce  ip  -  -  1 

in  the  equations  (6a),  (6b)  and  (8)  which  can  be  rewritten  in  the  following  manner 
2 

(for  slow  wave  0  >  e2  >  1): 

.2 


k  (1  -y)  f2 

u= - ip*(l -x)  +  2<py  +  y* 

x  +  y -  1  L  J 


y 

k2 


v  = 


-  <P  Ixy  -  2(x  +  y  -  1)1  +  2y(l  -  y) 

2(x  +  y  -  1)  L  l  J 


(12a) 


(12b) 


2 

v  -  u  = 


2 

k  y 


L  2(x  +  y  -  1) 


ip2  x2  +  Mfl/  +  1)  x(l 


-y) 


<p+  1  =  0  > 


(12c) 


(13) 


Case  I 


u  <  0 


rf2  >  0  ,  rf2  <  0 

1  2 


This  situation  can  be  satisfied  in  the  following  three  ways. 
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1  -x>  0 
1  -y>  0 
l>x  +  y-  l>  0 


(14a) 


(14b) 


1  -  x  <  0 

1  -  y  >  0  (14c) 

<P<  — y-r  (i+  fx) 

Although  the  condition  in  (14a)  does  not  explicitly  depend  on  ip,  for  slow  wave  ip 
must  satisfy  the  condition  (13).  On  the  other  hand  the  conditions  in  (14b)  and  (14c) 
show  that  besides  the  restriction  imposed  on  x  and  y,  ip  must  satisfy  two  simul¬ 
taneous  conditions.  More  precisely,  the  inequalitites  in  (14b)  suggest  that  ip  must 
be  greater  than  a  certain  value,  i.  e.  the  wave  propagation  which  satisfies  (14b)  is 
possible  for  a  value  of  ip  above  a  certain  value,  hi  other  words  this  passband  of 
wave  propagation  depends  on  the  degree  of  slowness  of  the  waves.  On  the  other 
hand,  the  condition  in  (14c)  can  be  met  only  for  a  definite  range  of  ip,  namely 

e2  -  1  <  ip  <  — - —  (1  +  Jic) 

x-l 
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where 


under 


For  v 


1  -  x  <  0,  and  1  -  y  >  0. 

Before  investigating  Case  II  and  Case  W,  it  is  desirable  to  find  conditions 
2  2 

which  u  >  0,  v  -  u  >  0,  v  >  0,  and  v  <  0,  respectively. 

For  u  >  0.  one  finds  the  following  possibilities 


1  -y<0 
1  -x>  0 


1  -y<0 
1  -x<0 


^  < 


x  —  1 


1  -  y  >  0 
1  -  x<  0 


0  > 


x  —  1 


> 


(1+  Vx)  J 


(1  +  fx)  J 


1  -  x>  0 

1  -  y  >  0  k 

x  +  y  -  1  <  0 


-  u  >  0,  the  following  conditions  must  be  satisfied 
1  -y>  0 


(15a) 


(15b) 


(15c) 


(15d) 


(16a) 
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(here  one  may  choose  either  x  +  y  -  1  >0,  orx  +  y-1  0) 


1  -y  <  0 


f  >  -11  +  2  ^(y  -  l)(x  +  y-  1) 


The  following  three  situations  satisfy  V  >  0: 


x  +  y- 1  >  0 


1  -y  >  0 


xy  <  2(x  +  y  -  1) 


U/  <  -  2^  ‘  -  g - 

I  2(x  +  y  -  1)  -  xy 


x  +  y  -  1  >  0 


i-y>  0 


xy>  2(x  +  y-  1) 


l  -y<  o 

xy  >  2(x  +  y  -  1) 


T  > 


xy-2(x+y-  1) 
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Finally  one  obtains  the  following  four  possibilities  for  V  <  0: 


1  -  y  <  0 


(18a) 


2(x  +  y  -  1)  >  xy 


1  -y  <0 

xy  >  2(x  +  y  -  1) 


xy  -  2(xfy-l) 


x  +  y -  1  >  0 

1  -y  >  0 

2(x  +  y-  1)  >  xy 

v  2(x  +  y  -  1)  -  xy 


x  +  y  -  1  <  0 
1  -y>  0 


(18c) 


(I8d) 


To  satisfy  the  requirements  for  Case  n  and  Case  in,  it  is  necessary  to 
satisfy  inequalities  (11b)  and  (11c)  respectively.  It  can  be  shown  by  a  little  analysis 
that  the  following  are  the  conditions  by  which  Case  n  and  Case  III  can  be  realized. 
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This  situation  can  be  obtained  if  the  following  conditions  are  met: 


1  -y  <  0 


1  -  x  <  0 


2(x  +  y  -  1)  +  x  {xy  -  2(x  +  y  -  1)J  >  xy  /x"  +  2(x  +  y  -  1) 


>  xy  +  2(x-  1)  •  i/(y-  1)  (x  +  y-  1) 


The  inequalities  in  (19)  are  obtained  from  (15b),  (16b)  and  (17c).  There  is 
no  other  possibility  which  can  satisfy  Case  II. 


In  this  case  one  can  show  that  there  are  only  four  possible  conditions  as 


1  -y  >0 
1  -  x  <  0 

4>  >  (1  +  sfx) 


(20a) 


Note  that  this  condition  is  the  same  as  (15c),  which  also  satisfies  conditions  (16a) 
and  (18c)  automatically. 
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1  -  x  >  0 

1  -  y  >  0  •  (20b) 

x  +  y  -  1  <0 

The  condition  (20b)  is  the  same  as  (15d)  which  also  satisfies  conditions  (16a) 
and  (18d). 

1  -  y  <  0 
1  -  x  >  0 

^  ~  \/(y  -  l)(x  +  y  -  1) 

This  condition  (20c)  is  obtained  by  combining  the  conditions  (15a),  (16b)  and 
(18a).  It  may  be  noted  that  condition  (15a)  automatically  satisfies  (18a)  also. 

The  fourth  possibility  of  realizing  Case  mis  obtained  by  combining  the 
conditions  (15b),  (16b)  and  (18b).  Since  it  is  not  obvious  that  these  three  conditions 
can  be  satisfied  simultaneously,  it  is  necessary  that  they  must  meet  the  following 
requirements  (a  detailed  analysis  is  omitted  for  the  sake  of  brevity) . 

1  -y  <  0 
1  -  x  <0 

xy  >  2(x  +  y  -  1)  r  (20d) 

2(x  +  y  -  1)  +  xy  \/x  >  xy  +  2(x  -  1)  \/(y  -  l)(x  +  y  -  1) 

>  2(x  +  y  -  1)  +  x  -jxy  -  2(x  +  y  -  1)J 
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Thus  it  is  found  that  when  the  condition  (10)  and  any  of  the  following  condi¬ 
tions  (14),  (19)  and  (29)  are  satisfied  simultaneously,  one  obtains  maximum  pass- 
band  of  slow  wave  propagation.  These  conditions  are  necessary  for  slow  waves. 

Study  of  Dispersion  Relation  for  Slow  Wave  Propagation  Under  Various  Special 
Situations 

Although  the  various  results  together  with  dispersion  relations  obtained  in 
the  preceding  chapter  are  valid  for  any  arbitrary  angular  variation  of  the  magnetic 
current  ring  source,  in  this  section  only  those  dispersion  relations  which  are  inde¬ 
pendent  of  angular  variation  (i.  e.,  n  -  0,  for  constant  amplitude  of  the  ring  source) 
will  be  considered. 

Since  the  solution  of  the  dispersion  relation,  appropriate  for  any  particular 
case,  together  with  the  expressions  in  (5)  for  rj'i2  and  rj'a2  gives  exact  information 
of  the  propagation  of  waves,  and  as  this  solution  cannot  be  obtained  analytically  in 
general,  the  information  obtained  here  without  actual  solutions  will  give  only  nec¬ 
essary  conditions  for  slow  wave  propagation.  In  general,  the  actual  solution  can  be 
obtained  only  by  numerical  computation. 

^Static  limit:  This  static  limit  is  a  good  approximation  in  the  following 
situations: 

1)  circumference  of  the  plasma  column  is  much  shorter  than  the  wavelength 
of  the  operating  frequency 

+Trivel  piece  in  his  work  discusses  this  problem  in  detail  and  his  method  of 
solving  this  problem  is  different.  Here  his  results  are  obtained  as  a  limiting  case. 
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2)  for  extremely  slow  wave,  i.e.,  >>  1. 

This  dispersion  relation  in  this  case  for  n  -  0,  can  be  obtained  from  the 
relation  (3)  of  Appendix  C,  and  it  reduces  to  the  following  form: 


V»’i  <Mn’ia)  \  <*a)  KJftfb) + Iq  bjK^a) 

e2rf  J  (rj'ia)  =  I  (^)K  (tfa)  -  I  (tfa)K  (*b) 
O  0  0  o  o 


(21) 


where 


n’ia  »  -  (ez/er)  *a 


r?  -a  - 


(22) 


The  above  relations  show  that  the  dispersion  relation  and  hence  fields  do  not 
depend  on  n'a  la  this  limit.  Moreover,  in  this  limit  the  magnetic  current  ring 
source  excites  only  E-type  mode.  But  in  a  general  anisotropic  medium  character¬ 
ized  by  e  of  the  form  shown  in  (1),  pure  E-type  and  H-type  modes  do  not  exist, 
i.e.,  they  are  coupled  to  each  other. 

From  the  properties  of  modified  Bessel's  functions  it  can  be  shown  that  the 

right-hand  side  of  (21)  is  positive  and  greater  than  unity  (it  approaches  unity  as 

►  oo).  If  e  /e  <0,  rj\  is  real,  a  solution  of  (21)  is  possible, 
z  r 

Since  e  /e  =  <  0 

z  r  x  +  y  -  1 

either 

x  >  1 

(23a) 

y  >  i 
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or 

x  <  1 

y  <.  i 

1  4  x  +  y  4  2 

If  e  /e  >0,  rj'ia  <  0,  and  rr'i  i®  purely  imaginary,  the  dispersion  relation  (21) 

2  1* 

bo  co  mo  s 

JTJT  Ij  (tfa  jcjt  J  k  <*a) KlXb) + IQ(^b)K^fa) 

r  e2  I  (tf  a  F7C )  I  <rfb)K(/a)  - I  (tfa)K(tfb) 
o  r  z  r  o  o  o  o 

A  solution  of  (24)  is  possible  if  e  <  0,  c  4  0  (since  e2  >  0).  As  it  is  also  known 

r  z 

that  Ix(z)  <  I&(z),  for  any  z  >  0,  it  is  necessary  that 


In  particular  if  e2  =  1  (i.  e. ,  when  the  plasma  column  is  surrounded  by  air),  the 
inequality  (25)  reduces  to  the  following 


Finally  the  inequality  (26)  can  be  shown  to  be  equivalent  to  the  following  two 
passbands  for  slow  waves  when  the  dielectric  surrounding  the  plasma  column  is  air: 
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y  -  1  >  0 

1  -  x  >  0  - 

x  +  y  >  2 


y-  1  <  0 


1  -  x  <0 
x  +  y  <  2 


> 


(27a) 


(27b) 


It  may  be  noted  here  also  that  when  the  plasma  column  is  surrounded  by  a 
dielectric  e3  >  1,  the  passbands  are  reduced  further. 

When  b  =  a,  i.  e. ,  when  the  plasma  completely  fills  the  waveguide,  it  can 
be  shown  that  the  corresponding  dispersion  relation  (in  the  static  limit)  reduces  to 
the  following 

Jo(rj'i  a)  =  0  (28) 

It  is  easy  to  show  that  only  real  values  of  i?'i  can  satisfy  the  above  equation 
(28).  Therefore,  in  this  case  also  it  is  necessary  that  <  0. 

A  study  of  the  relations  (21)  and  (22)  reveals  that  for  x  =  1,  y  =  1  or 
x  +  y  =  1,  the  dispersion  relation  (21)  does  not  possess  any  solution,  which  is 
equivalent  to  saying  that  these  points  represent  cut-off  for  the  slow  wave  propaga¬ 
tion. 

It  may  be  remarked  here  that  the  passbands  for  slow  wave  propagation  when 
the  plasma  completely  fills  the  waveguide,  give  maximum  range  for  the  case  when 
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the  plasma  column  partially  fills  the  waveguide  (provided  is  real)  in  the  static 
limit.  These  maximum  passbands  are  depicted  in  the  following  Figures  3  and  4, 
subject  to  the  conditions  (23a)  and  (23b). 


n'i 


+u*  -  u2)wa 

±/(u^  -u2)(w2  -ua) 
r  n  c 


(29)+ 


Since  group  velocity  is  defined  as  dw/dJf,  Figures  3  and  4  show  that  this 
value  can  also  assume  negative  values  —  which  proves  the  existence  of  backward 
wave  in  such  a  structure. 

With  a  few  more  remarks,  the  discussion  of  the  static  limit  case  will  be 
concluded.  The  static  limit  results  are  reasonably  valid  for  extremely  slow  wave 
propagation,  as  pointed  out  in  the  beginning  of  this  section.  In  this  limit  r)'j  does 
not  appear  in  the  dispersion  relation,  showing  that  the  field  components  for 
extremely  slow  wave  propagation  do  not  depend  on  r/j,  when  the  source  of  excita¬ 
tion  is  a  magnetic  ring  current.  In  other  words  in  this  limit  an  H-type  mode  is 
not  excited.  This  does  not  mean,  however,  that  in  this  limit  f] '  =  0.  In  fact, 

TJ '  =  -  j  X,  a  large  imaginary  number,  which  shows  that  a  wave  dependent  on  fj'2 
decays  away  very  rapidly.  So  it  may  be  conceived  that  in  this  limiting  condition 
the  H-type  mode  is  very  weakly  coupled  with  the  E-type  mode  and  the  components 

representing  H-type  mode  are  also  very  highly  attenuated.  It  may  be  noted  here 
“ - -  - - - 

+  In  8  ,  Trivelpiece  has  also  obtained  similar  results. 
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also  that  the  results  for  ^/ej  >>  1,  essentially  correspond  to  those  for  static  limit 


General  Dispersion  Relation  (when  n  =  0)  for  Slow  Wave  Prop»ff»Hnn 

It  can  be  shown  from  equation  (561)  of  Chapter  1,  that  for  slow  waves  (when 
H  =  _j  &  >  S  >0),  the  dispersion  relation  becomes 

ez^'a2  “  v\2)>f0(a)Go(a.)J1(r)\a.)J1(n'2&)  Cafa'a2  _  n'i2)Crt(a)S  (a)J  (rj'xa)J  (rj*2a) 

- - - -  +  - 2-  .  o - 2 _ 

k*(0*-  e2) 


n' i  n’a 


-^(rj'jaJJ^rj'ja)  -  _  _ 

[ea  M  Co(a)GQ(a)  -  *0(a)So(a)J 


kn. 


where 


J0(n'2a)Ji(r?'ia) 

+  TFT,  bs  Co(a)Go(a)  ezMjo(a)So(a)J  -  ® 


n2  «  -S*  =  k’ea-X5 

tf/k  =  |3 


M  -  k»(e2/er)(er  -  ?)  -  ^ 


S  = 


k*e 

~r  K  -  ^  -  n’i2 

r 


J  (a)  =  I  (<£  b)K  (f  a)  -  I  (S  a)K  ( S  b) 

O  0  O  0  0 


GQ(a)  =  I1(4b)Ko(tfa)  +  Io(6a)K1(5b) 


(30) 

Ola) 

(31b) 

(31c) 

Old) 

Ole) 
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S  (a)  *  Ii (6  b)K(^a)  -  sl)K.(£  b) 

o  1  A 

Claj  *  a)Ko(  &  b)  +  I  (S  b)K^  a) 


(31g) 

(31h) 


In  connection  with  the  solution  of  the  dispersion  relation  (30),  appropriate 
for  slow  waves,  no  general  discussion  can  be  made.  Only  a  numerical  solution 
subject  to  the  expressions  for  rj'i2  and  rj'8*  given  in  (5),  can  give  the  actual 
nature  of  slow  wave  propagation. 

The  dispersion  relation  (30)  will  be  solved  subsequently  for  a  special  case 
stated  in  (14a),  for  which  u  <  0,  t)'i  Is  real  and  n'j  i®  Imaginary. 

For  Zero  Magnetic  Field  (with  n  =  0) 

It  been  shown  in  Appendix  C  that  in  this  special  case  the  dispersion 
I  relation  ms  the  following  particular  form  (for  surface  waves) 


S  e  JiWia) 

z 


n'l  erMn'ia) 


■where 


I  (S  b)K  XI  a)  -  I  (£a)K  b) 

O  0  0  0 

(32) 

n't’  *  *  -P* 

(33a) 

-  i * 

(33b) 

|+This  result  agrees  with  that  obtained  in  ^7^  . 
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In  this  case  of  isotropic  plasma  only  an  E-type  mode  is  excited,  due  to  the 

type  of  the  source  of  excitation  chosen.  Since  e  is  always  less  than  1  and  >  1. 

z 

it  can  be  shown  from  (33a)  and  (33b)  that  when  rj  is  imaginary,  n'i  Is  also  imag¬ 
inary.  In  this  circumstance  a  surface  wave  is  possible  if  e  0,  moreover,  since 

2  ^(pa) 

the  left  hand  side  of  (32)  is  greater  than  unity,  p/S  <  1,  and  y  <  1,  it  is 


u 


also  necessary  that  /e2  >  1,  i.e.,  u  ^^5+^ 


For  Infinite  D.C.  Magnetic  Field  in  the  z-Dlrection 

In  this  case  the  dispersion  relation  for  slow  waves  has  the  same  form  as 

(32),  with  rj'i2  =  €zrtl  =  *«2  &  2  •  It  can  be  shown,  ,  that  slow  wave  is  possible 

if  e  <  0,  which  makes  r)'i  real.  For  zero  or  infinite  d.c.  magnetic  field  in  the 
z 

axial  direction,  an  elaborate  investigation  has  been  made  in  Q7J  . 


Propagation  of  Slow  Waves  in  an  Infinitely  Long  Column  of  Plasma  Embedded  in  an 
Unbounded  Medium,  with  an  Axial  Uniform  Static  Magnetic  Field 

To  investigate  all  related  results  a r  '  the  nature  of  slow  wave  propagation, 
in  this  case,  it  is  only  necessary  to  let  b  H  00  in  the  corresponding  results  of 
the  waveguide  problem,  with  r?  =  -82  =  k ae2-^Ca.  In  this  case  the  dispersion  rela¬ 
tion  takes  the  following  form  which  is  equivalent  to  equation  (22b)  of  Appendix  C 

(with  u  =  u  •  1,  p'  *  0,  H2  =  1) 
r  z 
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ez  c2  K2($a) 

17*1 17*2  J  <n'ia)J  (rj'aa)  +  $2  K2(£a) 

0  0  o 


e2  Kj($  a)  p  Ji(r)'ia)  Ji(r)'2a) 

nVa^W-n'i2)  p’2*5  Jo(rj’ia)  '  n'lM  JQ(<i;  a)  _ 

ez  K a)  p  Jitn'ja)  Ji(r)'ia) 

8  (r)'2S-r)’i2)  KjSa)  plS  JQ(rj’2a)  17  J^rj'ia) 


(34)+ 


The  following  identities  are  found  useful 

^(1  -  e  XSrj'i3  +  Mn'i2)  =  e,2e  (k2  +«*2)(k2e  +  **) 
r  z  z  z 

^(1  -  ezKSr7'22  -  Mr)*!2)  =  e^S  -  M)  -/*)(! +  er>  -  kV2 


(35a) 

(35b) 


To  derive  relation  (34)  it  also  has  been  assumed  that  ■  0,  which  may  be 
interpreted  as  following  from  taking  the  excitation  to  be  a  constant  nmgnetic 
current  ring  source. 

When  the  d.  c.  magnetic  field  is  either  zero  or  infinity,  only  the  E-type 
mode  is  excited  due  to  the  type  of  source  chosen  here.  In  this  particular  case  the 
dispersion  relation  takes  the  following  form 

+This  result  agrees  with  the  corresponding  result  in  ^4^  ,  when  the  identities 
(35a)  and  (35b)  are  used  (with  e2  =  1).  It  should  be  noted,  however,  that  to  obtain 
radiated  fields  in  this  configuration  the  present  limiting  process  is  not  valid. 
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^•Wia) 


Kj(5a) 

K  ( I  a) 
o 


If  the  axial  static  magnetic  field  is  zero,  €  -  e  ,  and  *j*i*  s  ~^a>  the 

r  z  z 

condition  of  slow  wave  is  exactly  the  same  as  stated  in  connection  with  the  similar 


situation  in  the  metallic  wave  guide,  namely,  e  <  0,  and  |e  |  /e2  >  1 . 

Z  I  z 

On  the  other  hand,  if  the  magnetic  field  is  infinity,  a  slow  wave  is  possible 

if  e  <0  and  r)\2  =  e  S2. 
z  z 

Static  limit:  It  can  be  shown  that  in  the  static  limit  with  b  -f  oo,  the  dis¬ 
persion  relation  (21)  reduces  to  the  following  expression 


Wi*)  *!<*») 


•Un'ia) 


KQ(/a) 


where 


n'i2  *  -  f  *a 
r 


If  c  /e  <0,  then  t)\  is  real  and  a  solution  to  (37)  is  possible  for  a  slow 

Z  1* 


wave. 


If  e  /c  >0,  then  17*1  is  imaginary,  i.e.,  17*1  =  -j P,  where  p  >0.  In  this 

Z  I* 

case  (37)  transforms  to  the  following  form 
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e^I^ap) 

Cjtfl^ap) 


Kj(^a) 

KjXa) 


Kj(^a)  Map) 

SinCe  k  0(a)  >  1  and  j  (ap)  ^  J»  14  iB  necessary  in  order  for  (39)  to 

have  any  solution  that  e  <  0  and  also  —  ^  )  1 .  It  has  already  been  stated 

r  ^2 

above  that  e  /e  >  0,  therefore,  e  <  0  and  e  <0. 
z  r  z  r 

If  one  writes  =  1-y,  as  defined  in  (7),  then  the  conditions 


|e  \P 

^0/  e  ^  0  and  rl  ^  1,  become  equivalent  to 


(y  -  l)(x  +  y  -  1)  .  2 

1  -x  ?  *« 

y- 1>  o 

1  -  x  >  o 

J  (40) 

Note:  The  discussion  on  page  497  of  QiJ  of  the  situation  where  u  =  0,  i.e. , 
r)'a  =  0  seems  to  be  inconsistent.  The  first  reason  is  that  when  rj'2  =  0,  it  can  be 
shown  from  the  general  expression  for  and  appearing  in  Chapter  I,  as  well 
as  in  Appendix  A,  that  electromagnetic  waves  which  can  exist  under  such  a  situation 
are  TEM  only.  In  this  statement  it  is  also  assumed  that  the  diagonal  components 
of  j£  are  finite  and  non-zero.  But  a  TEM  wave  cannot  exist  in  a  structure  con¬ 
sidered  by  the  authors  of  ^4^  .  A  similar  inconsistency  appears  on  pp  183  -  185 
of  [e]  ,  discussed  by  Agdur. 
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Secondly,  the  authors  of  QiJ  consider  a  case  when  y-+ao  (and  rj'a  =  0)/ 
and  simplify  the  dispersion  relation  to  an  expression  containing  a  logarithmic  term, 
although  the  modified  Bessel's  functions  K^(z)  and  K^(z)  do  not  behave  as  logarithmic 
functions  for  large  arguments. 

Any  other  interesting  situation  can  be  studied  by  considering  the  corres¬ 
ponding  dispersion  relation  given  in  Appendix  C. 

Expressions  for  E  Which  is  Independent  of  Angular  Variation  (i.  e. ,  n  =  0) 

_ _ z _ 

Since  E  plays  an  important  role  in  a  plasma,  its  expression  will  be  given 
z 

here  for  n  =  0.  |E  I  will  also  be  calculated  numerically  as  a  function  of  r  for  a 


special  case  of  slow  wave. 

$  $ 

-  .  •  ?21C,oW  f. 

Ez  =  -  2  2_>  .  LS 

,5liV"2ir)  -  MlV”’‘r)] 

X*1 

for 

0  4  r  4  a  (41) 

c  Si5liJ.(”21a)  '  Ml  W 

21  '"21  •  "'ll*  4o<‘> 

(42a) 

Si  ■  T  -"u 

r 

(42b) 

Mi  =  f 

(42c) 

r 
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Liv,  |”ii| 8  [v4  “A  - u  v'lM'] 
liVz^i5ii  K-’Vi'V'  •  RiViisiJ 

LiErtzn'li,,2i?li  fr’,2iM|Tl  +  wo*le'n'llili  J 

-  LlV,  |5lif  'rVllSl] 

[v> w*  r.<-,»r-o n* iwr] 

■  ^  [v>t?2lfW)fy’lb)  +  «o|53^lSb)Jl<"lb|] 

[*•*  ^rm-oN1  fivf] 

jwc  e  c'X.  r  “I 

,,2  ‘  ■ -  Lvo^  -  j„<v>J 

w21  -  ”  li)ErGo<a) 

AVX.e' 

_ o  z  i _ 

l'>lin2l<n2f  "lA 

erMi(kS'r  "*?  ‘ 


we  e  e' 
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,  e  M.X  -X4e  (l^e  -X*) 

_/  r  1  1  1  z  r  1 

Ri  =  e 

z 


/l0(r)  *  W’W  ■  WW 


SiQ(a)  =  J^WN^a)  -  Jj^Nj^b) 


(43d) 


Gio(a)  =  -  Jo(r]ia)N1(r)ib) 


Cio(c)  =  Ji(n1c)No(nib)  -  J/^WNifyc) 


\  (n’21r)rdr 


r)rdr 


Ji(n'2iDJi  (rj^rjrdr 


Ji(n’2r)Ji  (n'2ir)  rdr 


Xi5  =  \  Ji^ir>Ji^ir)  rdr 


Ijg  =  \  rdr 


(44b) 


(44c) 
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ni 

PROPAGATION  OF  SLOW  WAVES  IN  AN  ANISOTROPIC  FERRITE 

In  this  chapter  the  general  problc  m  solved  in  Chapter  I  will  be  applied  to 
the  study  of  wave  propagation  in  an  infinitely  long  anisotropic  ferrite  column 
enclosed  in  a  dielectric  medium,  which  is  again  enclosed  by  a  perfectly  conducting 
metallic  cylindrical  waveguide.  All  other  conditions  are  similar  to  those  described 
in  the  previous  chapters.  In  a  ferrite  medium  with  an  axial  static  magnetic  field, 
anisotropy  is  exhibited  by  the  following  dyadic  form  £2  J  of  n 


rr 


JM 


0r 


r0 


69 


zz 


(1) 


where  ^  -  nQQ  -  ^  =  1 


"re  =  "0r  =  =  - 


oP 

'  7^or 

p 

l-o* 


JLt  =1 

ZZ 


P  =  y 


M 
_ c 

H 


0 

h  ~ 


(2a) 

(2b) 

(2c) 

(3a) 

(3b) 
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=  d.  c.  magnetization 
Hq  s  d.  c.  magnetic  intensity 
7  =  gyromagnetic  ratio  for  electron 

It  can  be  shown  from  (2)  that  n  and  n'  are  related  by  the  following  relation 

r 

u  -  1  =  OM'  (4) 

r 

It  will  be  assumed  that  the  relative  dielectric  constant  of  the  ferrite  is  <  j , 
a  scalar  quantity  and  the  medium  surrounding  the  ferrite  has  relative  dielectric 
constant  C2  and  relative  permeability  1. 

With  the  above  assumptions,  the  magnetic  fields,  dispersion  relations  etc. 
for  this  problem  under  any  limiting  conditions  can  be  easily  derived  from  the  cor¬ 
responding  results  given  in  Chapter  I  and  Appendix  C .  Therefore,  no  detailed  dis¬ 
cussion  will  be  given  in  this  chapter. 

It  may  be  noted  here  that  when  a  ferrite  column  or  a  plasma  column  is 
situated  in  an  unbounded  dielectric  medium,  the  boundary  conditions  for  E  and  H 
in  both  cases  are  identical,  consequently  any  general  expression  for  one  situation 

can  be  derived  from  the  other,  using  the  duality,  provided  n  is  not  replaced  by 

z 

unity  in  any  general  expression. 
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Expressions  for  Transverse  Wave  Numbers 

If  the  values  for  £  and  €  given  above  are  substituted  in  equations  (34)  and 
(35)  of  Chapter  I,  one  obtains  the  following  expressions 


.2 

1 1,2 


V  +  \/va  -  u 


(5a)+ 


-  jjp  |k»C,(o(lr -dl-jAr} 


fr  ")  a  ""I1 

+  2^-  ^ka€1(oMr-#i,)-<a<Tj>  +  4ka*a«1  1 


1*1  + 


/2 


(5b) 


(6a) 


2u 


U  =  -J-  [oftiii  -Jf1)*- Afu4] 
’  •» 


(6b) 


For  slow  surface  waves  the  following  condition  should  be  satisfied 


—  •  P  >  «2 

k2 


1  *  -  J*  .  *  >  0 


(7) 


I  T  > 

+  For  ferrite  problem  the  relation  »}i  =  cannot  be  satisfied,  since  V  -  U>  0. 
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(2  t2 

Since  =  1,  equations  (5)  show  that  V3-  U  >  0  and  and  *)  2  are  real 
for  real  p,  6!  and  €  2.  Therefore,  no  instability  phenomena  appears  in  the  case 

I  t 

of  a  ferrite  column.  It  may  be  noted  that  j  and  *)2  may  assume  purely  imaginary 
values  as  follows: 


Case  1 


Case  n 


Case  in 


»2 

Vi 


»2 

^2 


»2 

^1 


i2 


if  U<0,  V>0,  or  V<0 


if  U>0,  and  V>0 


if  U>0,  and  V<0 


(8a) 


(8b) 


(8c) 


The  above  information  and  the  general  results  given  in  Chapter  I  and 
Appendix  C  are  sufficient  to  obtain  any  particular  result  for  a  ferrite  column. 

It  may  be  noted  that  an  electric  current  dipole  source  is  more  appropriate 
for  a  ferrite  problem  than  a  magnetic  current  ring  source.  The  field  expressions 
for  an  electric  dipole  source  can  be  obtained  easily  by  using  the  appropriate  dyadic 
Green's  functions  developed  in  Appendix  B. 
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IV  CONCLUSIONS 

In  conclusion  it  may  be  mentioned  that  the  work  described  here  gives  a 

I 

systematic  rigorous  approach  of  solving  a  source  problem  (not  necessarily  a  ring 
source)  involving  a  homogeneous  anisotropic  cylindrical  structure  bounded  by  con¬ 
ductors.  Since  the  source  free  solutions  are  capable  of  representing  all  possible 
modes  for  the  structure  of  the  problem,  total  fields  due  to  any  arbitrary  source  of 
any  kind  (namely  electric  or  magnetic  current  source)  can  be  calculated  by  using 
the  appropriate  dyadic  Green's  function.  If  there  is  more  than  one  source,  the  total 
fields  can  be  obtained  by  using  the  superposition  theorem,  provided  there  are  no 
interactions  among  the  Individual  sources.  As  the  present  analysis  considers  a 
magnetic  current  ring  source  of  arbitrary  angular  variation,  the  results  can  be  used 
for  any  given  angular  variation  of  the  source. 

From  the  general  dispersion  relation  which  is  an  eigenvalue  equation  and 
independent  of  source,  various  interesting  special  cases,  some  of  which  are 
already  known,  have  been  studied.  The  limiting  proct  ,*es  used  in  obtaining  the 
dispersion  relations  for  these  special  cases  can  also  be  used  to  obtain  the  expres¬ 
sions  for  the  fields  in  the  corresponding  situations. 

The  analysis  for  the  plasma  problem  which  is  a  special  case  of  a  general 
anisotropic  medium  characterized  by  dyadics  C  and  n,  emphasizes  the  slow  wave 
propagation.  Here  the  necessary  conditions  for  slow  wave  propagation,  including 
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a  number  of  special  cases,  have  been  obtained  from  the  expressions  for  the  trans¬ 
verse  wave  numbers.  These  necessary  conditions  also  give  maximum  passbands . 
The  sufficient  conditions  and  henoe  the  actual  passbands  can  be  obtained  from  the 
solution  of  the  dispersion  relation.  For  slow  wave  propagation,  some  passbands 
depend  on  the  degree  of  slowness  of  the  waves.  The  degree  of  slowness  of  the 
waves  depends  on  the  relative  dielectric  constant,  €  j,  of  the  medium  surrounding 
the  plasma  column,  when  all  other  parameters  are  kept  constant.  It  has  been 
shown  that  the  higher  the  value  of  €  2,  the  slower  the  phase  velocity  of  the  wave. 

In  other  words,  for  a  given  €  3  there  is  a  minimum  phase  velocity  for  which  a  cor¬ 
responding  slow-surface  wave  can  propagate.  Most  of  the  energy  of  the  slow  waves 
considered  here  is  confined  into  the  anisotropic  plasma.  In  general,  the  lower  the 
phase  velocity  of  the  surface  wave,  the  lower  the  amplitude.  Not  all  the  various 
passbands  for  slow  wave  propagation  in  an  anisotropic  plasma  column,  mentioned 
above,  are  known  in  the  literature,  at  least  to  the  best  knowledge  of  the  author. 
Although  these  passbands  could  be  obtained  without  any  consideration  of  the  pres¬ 
ence  of  a  source. 

In  the  case  of  an  unbounded  homogeneous  anisotropio  medium  where  a  TEM 
wave  can  propagate,  conditions  of  TEM  wave  propagation  in  the  direction  parallel 
to  or  perpendicular  to  the  static  magnetic  field  are  obtained  from  the  general  ex¬ 
pressions  of  the  transverse  wave  numbers.  It  is  also  shown  that  the  condition  of 
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TEM  wave  propagation  in  the  direction  of  the  static  magnetic  field  is  equivalent  to 
the  vanishing  condition  of  the  product  of  the  transverse  wave  numbers.  This  also 
establishes  the  fact  that  the  consideration  of  the  situation  under  which  the  product 
of  the  transverse  wave  numbers  vanishes  is  not  justified  in  connection  with  the 
wave  propagation  in  a  bounded  medium  which  cannot  support  TEM  waves.  In  other 
words,  if  a  bounded  isotropic  medium  cannot  support  TEM  waves,  so  also  is  the 
case  for  a  bounded  anisotropic  medium.  For  not  being  able  to  recognize  the  fact 
that  the  condition  of  TEM  wave  propagation  in  the  direction  of  the  static  magnetic 
field,  is  equivalent  to  the  zero-value  of  the  product  of  the  two  transverse  wave 
numbers,  some  authors*  discussed  the  possibility  of  wave  propagations  in  a 
bounded  anisotropic  plasma  column,  under  the  situation  for  which  the  product  of 
the  two  transverse  wave  numbers  vanishes. 


+  Agdur,  pages  183-185  of  Ref.  and  the  authors  of  Ref.  ^4^  ,  page  497. 


-  THE  UNIVERSITY  OF  MICHIGAN  - 

4386-1-T 

APPENDIX  A 

MAXWELL'S  EQUATIONS  FOR  ANISOTROPIC  MEDIUM 

The  medium  to  be  described  here  is  characterized  by  a  relative  dielectric 
(permittivity)  tensor  €  and  a  relative  permeability  tensor  n  having  the  following 
particular  form 


The  above  representations  show,  in  both  cases,  that  the  transverse  com¬ 
ponents  and  longitudinal  components  of  the  tensors  are  uncoupled,  where  £  ^3  and 
I Ugg  correspond  to  the  longitudinal,  a  preferred  direction  (say  z-direction)  compo¬ 
nents  of  the  tensors  £  and  u  respectively.  Therefore,  £,  and  n  can  be  written 
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In  the  following  manner  also 


£  *  £tt+VoC33 

(3) 

t  ,,‘tt+2oV33 

(4) 

where  and  are  tensors  (transverse  to  z-direction)  and  z 

tt  tt  — 0 

in  the  z-direction. 

is  a  unit  vector 

The  Maxwell's  equations  for  anisotropic  media  with  sources  have  the  follow¬ 
ing  form  (the  time  dependence  being  e^1*). 

V  x  E(r)  *  -  ju^(r)  •  H(r)  - 

(5) 

V  x  H(r)  *  jug  g(r)  *  E(r)  + 1  (r") 

*  o  —  ~e 

(6) 

V  ■  /i( r)  *  H(r)  =  0,  for  r  /  r' 

A/ 

(7) 

V  *  €(r)  ‘  E^(r)  =  0,  for  r  /  r" 

where, 

r  =  observation  position  vector  (3-dimensional) 

I  (r')  =  magnetic  current  source  at  r  -  r' 

“TO  ““ 

X  (r")  =  electric  current  source  at  r  =  r" 

(8) 
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In  the  following  discussion  the  transverse  vector  and  transverse  operators  (which 
are  designated  by  the  subscript  t)  correspond  to  any  plane  transverse  to  the  z- 
direction,  and  the  transverse  plane  may  have  any  arbitrary  cross  section.  That 
is,  the  following  derivations  are  suitable  for  any  cylindrical  geometry  having  the 
z-dlrection  as  its  axis. 

First  of  all,  it  will  be  shown  (see  (jLoJ  ,  [llj  ,  [12J  ,  and  [14J  )  that 

the  longitudinal  fields,  E  and  H  ,  can  be  expressed  in  terms  of  the  transverse 

z  z 

fields,  E.  and  H. .  Secondly,  it  will  be  demonstrated  that  the  transverse  fields 
t  t 

can  also  be  expressed  from  the  knowledge  of  the  longitudinal  fields.  In  the  parti¬ 
cular  problem  discussed  in  the  text,  the  latter  method  has  been  adopted  for  the 


solutions  of  Maxwell's  equations. 

If  equations  (5)  and  (6)  are  multiplied  by  zq  in  a  scalar  product  fashion, 
Ez  and  can  be  expressed  in  the  following  way 


Ej  *  >“‘0‘33 


X  z  - 


J“‘o*33 


H 

z 


J“V33 


I 

mz 


)UV33 


(10) 


where  I  and  I  are  the  z-components  of  I  and  I  respectively, 
ez  mz  — e  ~m 

Now  taking  the  vector  product  of  (5)  and  (6)  with  z  ,  one  obtains 
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z  x  ( V  x  E)  =  -  j  uji  z  x  n  '  H  + 1  xz 
“o  —  <r~o  “  -m  “o 


z  x  (V  x  H)  =  jwe  z  x€*E  +  z  xl 
— o  v  —  o-o  —  ~o  ~i 


Introducing 


V*  l  ■ 


and 

one  obtains 


E  =  E  +  z  E 
-  — t  — o  z 


H  =  H+z  H  , 
-  — t  -o  z 


and 


^  x  (V*  E)  =  VtEz-  £  ^ 


A>  *  <V*  ®  *  V.Hr  -  £ 


Now  equations  (11)  and  (12)  can  be  rewritten  as 


V.E  -  -f-  E* 
t  z  9z  — t 


-  Jwm  z  XJU.  H.+I.XZ 

o-o  M  -t  — mt  -c 


and 


V.H  -  -r—  R  =  jwe  z  x  6  ‘  E+ z  x  I 
t  z  3z  ~t  o-o  ~t  -t  -o 


— et 


(11) 

(12) 


(13) 


(14) 


where  1  and  j[0t  are  transverse  components  of  1^  and  I  respectively. 

8 

Again  operating  (14)  by  jwji  z  x  u.  and  (13)  by  —  from  the  left  one 

0“0  az 

obtains  the  following  expressions 
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to"  VtV  i?  = 


-  ju>M  T"(z  x  u  '  IT)  +  -7-  I  *xz 
o  9z  -o  ~t  — t  9z  -mt  -c 


(15) 


and  Jw#!^  x  <Ht  '  V  tHz)  -  jWM0  Tz  K  X  Hi  '  ^ 

■  -  x[«t '  *o  x  <*t '  v]  +  Ju Vo  x  (ift  '  — o  *  V 


(16) 


where  k2  =  u  e  . 

o  o 


Adding  (15)  and  (16),  one  can  show  that 


7  V  E  +  Ju#i  z  x  (u.  *  VH  )  +  —  z  x  I  *  +  (#i, ,  1*  ~  j  «£  x  1.)  *  I  * 
dz  t  z  J  o~o  *H  t  z  dz  -o  -mt  ll*t  12-o  ~t  -et 


■  *A-Jk,a3io  * 


where  1.  =  transverse  unit  dyadic 


a3  =  U1 1 *12  + ^12*11 


d3 


a4  *  ka<*nMu+  *«',U)+  3z 


<wnit  -  ^lA  a  V  a  <*t '  a  V 


(17) 

(18a) 

(18b) 

(18c) 

(18d) 


Taking  the  vector  product  of  (17)  with  zq,  one  obtains  the  following  independent 


equation: 
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H  -o  x  VtEz  -  '  V«Bz  '  Tz  *mt *  ““lA  +MlA>  1  -V  '  -et 


”  Vo  *  5t  +  )k>a3^t 


Eliminating  z  x  E  from  (17)  and  (19)  by  multiplying  (17)  by  a.  from  the 

left  and  (19)  by  j  k^a^  and  then  adding  the  results,  E^  can  be  expressed  in  terms  of 

E  ,  H  ,  I  .  and  I  in  the  following  way: 
z  z  ~ mt  — et 

pl^t  =  [-a4  A  VtEz  ‘  “Woa2  VtHzl  " -o  xDk’*3  #zVtEz  +  iU"oaiVtHz] 

-  [aUt  +  X  ij  '  -et 


-fa.  T"z  x  1.  -jk*a~  1*1  '  I 
L  4  9z  — o  -~t  3  9z  ~tJ  — i 


where  the  following  relations  have  been  used. 


,4  2  2 

pl  =  h  a3  -  a4 


ai  =  Vll  •  k2a3,,12  '  k2*U(*‘l2l  '  “l2)+l‘u  ~2 

oz 

2 

a2  =  k  V*ll”Vl2  =  k  612^M11  “  #<12>  “  M12  ^2 


(21b) 


and  the  identities 
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-oX  Ht  VtH,  =  "lA  x  V.  +  ,»*12\H. 


(22b) 


By  a  similar  method  (or  by  duality  of  (20) ),  one  can  obtain  the  expressions  for 
in  terms  of  H  ,  E  ,  I  .  and  I  .  in  the  following  form: 

Z  Z  nit  ™ ’ "0t 


p,H.  =  f-a.  —  V  H  +W6  a  V  E  1 
1-t  L  4  8z  tz  o  2  tzJ 

-z  x  Tjk^a,  “■  ^.H  -  jue  a,V*E  1 
-o  L1’  3  dz  tz  ol  tzJ 


-C*A+JVo  2iJ  lmt 

+  [a  ~r~  z  x  1  -  jk2a_  1.1 
L  4  9z  — o  ~t  3  9z  ~t-J 


where 


al  =  k  Mll(eil  "612)  +  €11  8z2 


.2  /-2  2  .  _  8 

a2  ‘  k  M12(€ll"  612)_  ei2  2 

dz 


(24b) 


Although  the  medium  is  anisotropic,  if  it  is  homogeneous  (i.e.,  components 
of  fe  and  fi  are  not  functions  of  position,  although  they  may  be  piecewise  con- 
stants)  and  source  free  (i.e.,  I  =  0  =  1),  one  can  show  that  £_  and 
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satisfy  the  following  equations  fusing  (20)  and  (23)  In  (9)  and  (10)  J  , 


'33 


,6  .  ^  x% 


£  +  -  OT  ^  s  ~  -  ~  o 

*  2  1  1  en“ii  3  * 


+  33 

2  «U<*U 


a 


%  . . 


1  z 


)uy33* 

eul‘ll 


a, 


3  z 


(25) 


(26) 


where  £  and  fr  are  solutions  of  homogeneous  (source- free)  Maxwell's 

Z  Z 


equations. 


To  obtain  the  above  two  expressions,  it  has  been  assumed  that  ~  = 
where  ^  is  the  propagation  wave  number  in  the  z-dlrection.  This  assumption  is 
permissible  in  the  situations  where  both  the  6  ,  j*  and  the  geometry  of  the  prob¬ 
lem  are  Independent  of  z,  subject  to  another  restriction,  that  tne  transverse 
anisotropy  of  j&  and  jU  are  not  coupled  to  the  longitudinal  anisotropy  of  €,  and  p 
respectively. 

To  obtain  solutions  for  £  and  ,  it  is  possible  to  have  4^  degree 

Z  Z 


equations 


in  £  and  ^ 


and  Jr  from  (25)  and  (26)  by  elimination.  Since  it  is  a  tedious 
z  z 


task  to  solve  such  equations,  one  can  alternatively  find  a  function  0  which  is  a 
linear  combination  of  <5  and  satisfying  a  two  dimensional  wave  equation. 


Let  such  a  choice  be 


t- 


(27) 
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Now  multiplying  (26)  by  j*  and  then  adding  the  result  to  (25),  one  obtains  the 
following  relations 

**»<£, +*fy+  —j-  [*[  -  “«0*v*] 

•nwii  L1  “  J  * 


(28) 


The  above  equation  can  be  represented  as  a  two  dimensional  wave  equation  in  0 
of  the  following  form 


7*  0  +  T2^  *  o 


(29) 


where 


*33  r.  .  1  2  **33  r  “w  1 

‘ii'-u  L*i'ue^J ' 1  L*»’  «  J 


(30) 


Solving  equation  (30)  for  a,  one  obtains 


®i,a  B 


*33*1  "  *zz*l  *  0*1*33 ~ai **33^  *  4k  ^*3*33  **33] 


1/2 


2uee  Jfa, 
o  z  3 


(31) 


Therefore,  the  roots  '  can  also  be  expressed  in  the  following  form: 

J2  C33*l  *33*1  ~  **33*1  .  0*1*33 iVW  ***  ^*3*33^33]  ^ 

,2  '  *11^1  " 


26  u 
11M11 


2*11/<11 


(32) 


The  equation  (32)  can  also  be  rewritten  in  the  following  form 


.2 


=  Vi/v2  -  U 


(33) 
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V 


1«33  +  V83 
2€11M11 


(34a) 


U 


£33^33 


€33/<33 
€1 1^1 1 


(34b) 


Relatione 

Between  Various  Parameters  Introduced  in  the  Above  Analysis 

First  of  all,  the  following  new  parameters  are  introduced  and  defined: 

S 

=  €33*i  .  y2  = 
'*11*11  1 

Mneii 

(35a) 

u 

.  *33*1  .v.,2  . 

W6o£33^a3a2 

XVI 

u  e  <2 

MU  11 

^11£11 

R 

*  */ff2  '  “Vi 

(36a) 

T 

■  “Vi  -  **4*1 

(36b) 

R» 

=  “*„*2“2  -  **4 

(37a) 

T» 

*  “*o*2tfl  '  **4 

(37b) 
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It  will  be  found  convenient  to  carry  out  algebraic  operations  in  the  sequel 
using  the  relations  listed  in  the  following  table.  In  tabulating  the  following  results, 
no  attempt  has  been  made  to  write  them  in  such  a  way  that  any  result  is  a  conse¬ 
quence  of  those  preceding  it. 


(38-1) 


*2.  =  €33al 
1  M11C11 


uy3jlVi 

'*11*11 


(38-2) 


£  a* 


U£  6, 


2  =  33  1  _  “o  33T“3“2 


121**3 


Mll€ll 


/ill€ll 


(38-3) 


*1  ”  °2 


ueo*l3S33 


(38-4) 


V*3 

or  or  =  -  1 

1  2  «.633 


(38-5) 


£  a* 

,  -SSL.  *.2 

^11^11  1 


ut6*3&al 

Mll6u 


(38-6) 


^  '  ^oa3^33 


(38-7) 
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Maj33!l_r2a 
pll*ll  2 


Mll€ll 


2  WW33 


a„  = 


al  "  °2 


al-°2 


S 


,2  ,2 

U  '71 


M 


r?  - 


alC33  *  ai  ^33 


*)!2  +  t);2  =  — — — — — 

^  /2  Vil 


S-  M  = 


MS  =  - 


p  =  k4a2  -  a2  =  -  €llMl1  t)  .2  t)|2 
1  ^  4  *33*33  ’*  ’i 


SR  +  MT  ft2) 


1-Me”'‘"  -  k 


R  =  3f»  «  -  up  al  «  — — —  — — 

4  2  0  2  “«o*3‘33 


33 


(38-8) 

(38-9) 

(38-10) 

(38—11) 

(38—12) 

(38-13) 

(38-14) 

(38-15) 

(38-16) 

(38-17) 


100 


THE  UNIVERSITY  OF  MICHIGAN 

4386-1-T 


T  =  u  V2  '  *Vl 


R'  =  USa2a2"^a4 


T'  =  u£0a2al-*a4 


k  fe33a3a2  ~  a4Seil^ll 
UV3633 


Ma2fell^U  a3a4633 


3* 


33 


Sa2  £1 1**1 1  "  *  a3a4e33 
^a3fe33 


633a:  ’  6 


1  "  61 1^11^2  =  £11MI1M 


€33al  “  eilMllT)l  =  ^l^ll8 

2  k2^*8e3a^3 


M33ai  "  ^11^11^2 


,2-2  2 
*  ^3^33^33 


-  ^11*11 


-SC  u 
11M11 


**  a3  V*ll  "  aiM12 


al  =  a4eil  '  k  a3fe12 


**  a3  a26ll  "  al612 


ai  =  V11  -  k  V12 


(38-18) 

(38-19) 

(38-20) 

(38-21) 

(38-22) 

(38-23) 

(38-24) 

(38-25) 

(38-26) 

(38-27) 

(38-28) 
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a2  ■  k  Yu  -  Vl2 

(38-29) 

*2  '  k  a3*U  "  a4*12 

(38-30) 

“Vi  '  ala4  =  "llpl 

(38-31) 

*  a4a3  '  *2al  "  *Vll*12 

(38-32) 

k%S  -  Vi  =eupi 

(38-33) 

“^•s  -  *lai  *  Vu‘u 

(38-34) 

ai  "  k  pll*ftll  ”a12>  ■‘’^*11 

(38-35) 

,  ,  /  2  2  Y  2 

ai  =  k  eilKl“Ml2>’XMll 

(38-36) 

a2  =  k  M12(eil  ‘  fe12)  +  ^612 

(38-37) 

a2  =  k  fe12^11  _M12)+;(2#i12 

(38-38) 

a3  °  ^llp12  +p11€12 

(38-39) 

a4  *  k2<pn6n+pi2si2)'‘<2 

(38-40) 
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It  may  be  noted  that  all  of  the  analyses  and  results  obtained  in  this  Appendix 
are  based  on  the  assumption  that  the  problem  has  a  cylindrical  geometry  with  axis 
in  the  z-direction  and  having  an  arbitrary  cross  section  which  is  independent  of  the 
coordinate  z. 

Propagation  of  TEM  Waves  in  an  Unbounded  Homogeneous  Anisotropic  Medium 

Since  the  foregoing  analysis  does  not  include  any  particular  boundary,  it  is 
valid  for  an  unbounded  medium  also.  Hence,  it  is  possible  to  obtain  conditions  for 
TEM  wave  propagation  in  a  direction  parallel  or  perpendicular  to  the  z-axis.  It 
may  be  mentioned  here  that  a  plasma  and  a  ferrite  with  a  static  uniform  magnetic 
field  in  the  z-direction  will  have  tensor  permittivity  and  tensor  permeability 
respectively.  The  forms  of  these  tensors  are  given  in  equations  (1)  and  (2).  As 
mentioned  earlier,  here  also  the  medium  considered  will  have  both  €  and  n  as 
tensors  with  constant  elements. 

TEM  Wave  Parallel  to  the  Magnetic  Field 

For  a  TEM  wave  in  the  z-direction  both  E  =0  and  H  =  0.  If  the  source 

z  z 

terms  in  equations  (20)  and  (23)  are  equated  to  zero,  non-vanishing  values  of  E^ 
and  H  are  possible  if  and  only  if  p  =  0,  when  E  =  0  =  H  ,  provided  the  elements 
e  jj,  £3g,  and  *igg  are  finite  and  non-zero.  The  condition  p^=  0,  gives  two 
TEM  waves  propagating  in  the  z-directlon.  These  two  waves  are  characterized  by 
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the  following  expression  of  the  propagation  wave  number  in  the  z-direction 

'  <6ii±6i2><"n±'‘i2)'  (39> 

k 

2  2  2  2 

Since  p  is  proportional  to  r| j  fj  ^  ,  the  condition  *)[  1 2  =  0,  is  equivalent  to 
TEM  wave  propagation  in  the  z-direction,  provided  the  diagonal  elements  of  € 
and  u  are  finite  and  non-zero. 

TEM  Wave  in  the  Direction  Perpendicular  to  the  Static  Magnetic  Field 

The  conditions  for  a  TEM  wave  propagating  in  the  direction  perpendicular 
(i.  e. ,  perpendicular  to  the  z-axis)  to  the  static  magnetic  field  can  be  obtained  upon 
substitution  of  X-  0  (i.e.,  ~~  =  0)  in  the  expression  (32).  This  substitution  gives 

<7Z 

» 

two  propagation  wave  numbers  -fj  ^  ^  which  represent  two  TEM  waves  in  the  trans¬ 
verse  plane  of  the  z-axis 
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|^1  J  by  an  entirely  different  approaoh. 

When  Pj=  0,  a  study  of  the  expressions  (33)  and  (34)  shows  that  *7^  =  0, 
provided  the  diagonal  components  of  €,  and  ^  are  finite  and  non- zero.  It  can 
be  shown  also  that  such  TEM  modes  in  an  unbounded  homogeneous  anisotropic 
medium  do  not  vary  in  the  transverse  plane,  but  in  a  coaxial  waveguide  TEM  waves 
behave  as  l/r  in  the  transverse  plane.  The  above  statement  follows  from  the  fact 
that  for  a  TE  M  wave  the  transversely  varying  part  of  the  transverse  fields  can  be 
derived  from  -  V^(g),  where  0Q>)  is  a  scalar  potential  dependent  on  the  transverse 
coordinate  £. 

The  two  waves  given  by  (39)  are  known  as  ordinary  and  extraordinary  waves 
in  the  literatures  of  the  ionospheric  wave  propagation.  The  permeability  of  the 
ionosphere  is  a  scalar  quantity  and  equal  to  that  of  free  space.  These  two  equations 
also  explain  the  phenomena  known  as  Faraday  Rotation.  The  two  waves  represented 
by  (40a)  and  (40b)  can  be  said  to  explain  £l]  the  phenomena  known  as  magnetic 
and  electric  Cotton- Mouten  effect. 
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APPENDIX  B 

CONSTRUCTION  OF  DYADIC  GREEN'S  FUNCTIONS 

* 

Although  a  construction  of  Dyadic  Green's  functions  from  the  source-free 
solutions  of  Maxwell's  equations  for  inhomogeneous  anisotropic  non-dlssipative 
media  in  a  uniform  waveguide  of  arbitrary  cross  section  bounded  by  a  perfect  con¬ 
ductor,  has  been  discussed  in  jjL  5tJ  ,  they  will  be  also  briefly  presented  here  for 
the  sake  of  completeness  of  this  work.  In  this  appendix  the  corresponding  results 
for  anisotropic  dissipative  medium  will  also  be  obtained. 

The  most  important  technique  involved  in  the  construction  of  dyadic  Green's 
functions  is  the  determination  of  an  appropriate  orthogonality  condition  among  the 
source- free  solutions  (i.  e. ,  eigenfunctions)  of  Maxwell 's  equations.  Methods  of 
finding  such  orthogonality  conditions  have  been  discussed  elaborately  by  the  authors 
in  ,  under  different  situations . 

Here  an  indirect  method  will  be  presented  for  the  construction  of  dyadic 
Green's  functions  [12J  .  In  this  method  it  will  be  assumed  that  the  sources  are 
due  to  some  discontinuities,  which  causes  discontinuities  in  the  fields  also. 

Dyadic  Green's  functions  Zfr,  r'),  T  (r,  r'),  Y(r,  r')  and  T  (r,  r') 

^  —  — em  “  ^  _  ~m©  — 

are  defined  by  the  following  expressions: 


5(r>  ■  -  Z<£.  r'l  •  I((rW  -  jjj  ^.1  '  IJrW  (1) 
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H(r)  * 


r') 


I  (r')dV'  - 
— m“ 


Tme(r,  r')  ■  Ie(r')dV' 


Where  the  E(r),  H(r),  I  ,  I  have  the  same  significance  as  given  by  the 
—  “■©  — m 

Maxwell's  equations  (5)  to  (8)  of  Appendix  A.  Instead  of  volume  currents,  if  X 

© 

and  I  represent  surface  currents,  the  volume  integrals  in  (1)  and  (2)  should  be 

"TO 

replaced  by  surface  integrals  (over  the  regions  of  surface  currents). 

In  the  following  are  given  the  physical  meanings  of  dyadic  Green's  functions: 


-  Z(r,  r')  ‘  u  =  electric  field  at  r  due  to  a  point  electric  current  source  at 

r',  directed  along  the  unit  vector  u . 

-  T  (r,  r')  '  v  =  electric  field  at  r  due  to  a  point  magnetic  current 

source  at  r',  directed  along  the  unit  vector  v. 


-  Y(r,  r')  *  v  *  magnetic  field  at  r  due  to  a  point  magnetic  current  source 
at  r'  directed  along  the  unit  vector  v. 


-  T  (r,  r')  '  u  =  magnetic  field  at  r  due  to  a  point  electric  current  source 
/^me  at  r'  directed  along  the  unit  vector  u. 


In  the  above  statements  the  point  source  means  a  source  which  has  spatial  variation 
as  a  Dirac  delta  function  8  (t  -  r'). 

Let  &  (r)  and  ^(r)  be  the  solutions  of  the  homogeneous  (source-free) 
Maxwell's  equations  (3) 

V  x  £  (r)  =  -  ‘  r)  (3a) 
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V  x  3^(r)  =  ju>  €q€(£)  ’  £  (r)  (3b) 

where  nig)  and  e  <£>  are  functions  of  transverse  coordinate  £  only. 

Since  under  appropriate  boundary  conditions,  £  (r)  and  ^(r)  form  a 
complete  orthogonal  set,  the  total  fields  E(r)  and  H(r)  due  to  any  arbitrary  source 
can  be  expressed  as  a  superposition  of  £  (r)  and  $(r)  in  the  following  way: 

E(r)  ■  y  '  Aa  l  Jr)  <4a) 


Hte)  ■  A  *J 


where  is  the  coefficient  of  expansion  corresponding  to  o- th-mode  (eigenvalue). 


Reciprocity  Relations  for  Homogeneous  Maxwell's  Equations 


To  establish  Lorentz's  reciprocity  relation  and  hence  an  orthogonality 
condition  it  is  desirable  to  consider  another  set  of  Maxwell's  equations.  This  new 
set  of  equations  is  sometime  s  called  the  Adjoint- Maxwell's  equations  .  After 
taking  complex-conjugates  of  these  so-called  adjoint  equations,  the  resulting 
Maxwell's  equations  have  the  following  forms: 


_  >•  +* 
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\7  *  *  £fl 


where  p+,  €+  =  adjoint  of  p,  and  respectively 


3 


(5b) 


complex  conjugate  of  the  transpose  of  ju  and  €,  respectively. 


.  .  U  =  P 

A/  A-' 


transpose  of  p 

**11 

~iM12 

*12 

M22 

0 

0 

transpose  of  £ 

611 

-j612 

j612 

6  22 

0 

0 

33 


(6) 


(7) 


331 


+,t  h 

Although  the  authors  in  [~9 1  have  used  the  symbols  £  R  for  £  R 

±  J  $  $  P  P 

i|4  j|/  I!  ^ 

and  „  for  Ar  in  general  c  Q  and  %  n  have  no  simple  relations  with 

P  P  P  P 

the  solutions  of  equations  (3a)  and  (3b).  Thus,  to  avoid  confusion,  the  symbols 
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and  have  been  used  here.  Here  &  Q  and  %'a  only  mean  the  solu- 

P  ”  P  P  P 

tlons  of  (5a)  and  (5b),  subject  to  some  appropriate  boundary  conditions. 

Now  multiplying  (3a)  by  and  (3b)  by  and  (5a)  by  and  (5b)  by 


(£  in  a  scalar  product  fashion  from  the  left  and  then  subtracting,  one  can  show 

*  IX  -  #8-  4  X  -b]'  ° 


—  Gt 

that 


Since 


and 


4 


(8) 


(9a) 


(9b) 


Let  U  =  7,  +  z  ~ —  ,  where  z  is  the  unit  vector  in  the  z-direction. 

Y  t  — o  8z  “o 


Now  using 


and 


■  v;.  i*; 


equation  (8)  can  be  rewritten  in  the  following  way 


■  La  *  *  0 = k '  a  X  a  *  ^  a] 


(10) 
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Now  integrating  over  the  cross-section  of  the  waveguide  one  obtains  (using  the  two 
dimensional  divergence  theorem) 

§  [£X  **  '  •2i]d8  = 


where  v  is  a  unit  outward  normal  vector  on  the  boundary  curve  s  of  the  waveguide 
cross  section  S. 

The  left-hand  side  of  the  above  expression  vanishes  on  the  boundary  of  a 
perfect  conductor 


.  the  orthogonality  relation  becomes 

nM;-*.-*.!.-*;] 


1)8  *  “a  . 


where  N  is  a  normalization  constant  and 
a 


&  a  =  1  ,  for  X  = 
ar(3  a  p 


=  0,  for  &a  f 


When  the  waveguide  has  a  reflection  symmetry,  i.  e. ,  when  the  properties  of  the 
waveguide  are  independent  of  the  coordinate  z,  the  orthogonality  relation  can  be 
rewritten  as 
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]]  *  Vs  *  nA3  ■  j&'A,  x  — B  118 


(12) 


or 


55  ^ ■  -to *  -o ds 


■  NA(S  *  Q  dS  (13> 


where  the  suffix  t  represents  the  transverse  components  of  the  fields. 

Construction  of  Dyadic  Green's  Functions  (see  jjlj  and  £12]  ) 

Knowing  the  total  fields  E(r')  and  H(r')  which  one  may  consider  are  due  to 

discontinuity  at  some  cross  section  S  of  the  waveguide,  one  can  write,  using 

% 


equation  (4) 


3^<E’>  *  ^  •  Sfe’>  ■  fy')  *  V  Z  A* 


0 


and 


1Q  *  *  S(r')  =  zq  x  £g(r’>  *  £ 


-O  ~fr-  '  >  a" or- 

a 


Integrating  the  sum  of  the  above  two  equations  over  the  cross  section  at 

and  using  (11),  one  obtains  the  value  of  the  coefficient  A  as 

a 


jj  [s(£’>  •#> 


)  x  z  +  H(r')  ’  z  x 

"O  - ““0 


dS' 


A  = 
a 


(14) 


2N 
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Therefore,  the  total  electric  field  at  any  point  r  due  to  E(r')  and  H(r')  Is  given  by 

S  ds' (l<£’>  '!£<£'>  xio+H(r')  •  zox$V]ga(r) 

*  E  '  E  - - 


2N 


which  can  also  be  written 


z 


z  x  H(r')  dS' 
o 


E(r')  x  z  dS* 
-  ~ o 


(15) 


Since  z^  x  H(r')  and  E(r')  x  zq  represent  sources  due  to  discontinuities 
in  fields  at  r',  one  can  represent 


z  x  H(r')  =  I  \r')  and  E(r*)  x  z  =  I  ,(r') 

“O  _  _et  -  — ^  — o  — mt  - 


Alternatively  one  can  consider  the  discontinuities  in  z  x  H(r')  and  E(r')  x  z 

— o  "  “  “O 

at  r*  are  due  to  actual  sources  let(£')  and  I  (rO  respectively.  Now  with  the 
above  assumption  if  one  compares  equation  (15)  with  equation  (1),  one  finds 


Z  (r,  r') 


y  *>> 

2N 

a  a 


(16a) 
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T  <r.  ,•>  .  £ 

A/em  *  ? 

a 


2N 


a 


Similarly  if  one  expresses  the  total  magnetic  field  as  H(r)  =  [  A  ^  (r)  and 

or 

uses  equation  (14),  following  the  above  procedure,  it  can  be  shown  that 

&<£>*>> 


r')  =  £ 


2N 


(17a) 


_  ,  ...  v 

2mei  -  ”  E  2N 

a  a 


(17b) 


Special  Cases 

It  has  been  pointed  out  before  that  in  general  there  are  no  simple  relations 

between  and  and  X  and  ^  ,  but  in  some  special  cases  these 
a  a  a 

relations  simplify.  For  example,  in  non-dissipative  anisotropic  media  «+  *  €, 
and  /u+  =  n  ,  i.e.,  6  and  n  are  hermitean  (self-adjoint)  dyadics.  In  this  case 

^  S  ^  ,  I*  *  ^  is  real. 

"or  a  a  a  a 


Therefore,  for  non-dissipative  medium,  the  dyadic  Green's  functions  can  be 
expressed  in  the  following  way 


Z(r,  r') 


Z 


*.<£>  K  <£’> 


2N 


or 


(18a) 
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(18b) 

(18c) 

(18d) 

Another  kind  of  orthogonality  relation  and  hence  dyadic  Green's  function  can 
be  constructed  in  the  following  way  (see  ).  These  results  are  particularly 
suitable  for  dissipative  medium. 

In  this  method  the  following  replacement  is  made 


H*  -¥  n  -  transpose  of  n  (19a) 

6  =  transpose  of  6^  (19b) 

$ 

#+  =  X  -  X  (19c) 

a  a  a 

it 

6“*  +  £'  (19d) 

a  a 

<19e) 


Due  to  the  transformations  given  in  (19)  and  the  Maxwell's  equations  (3)  and 
(5),  the  following  simple  relations  can  be  established 


T  (r,  r')  =  T* 
-'em'-  —  ‘ 

a 


6aV*J& 


2N 


Xfcr')  -  2-.  — i5 - 


‘  Z,  - 2N - 


me 


a 
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*12'  “12-  ■£>  *  +  -  ei2--“l2-  (20>) 

and 

612'  M12*  =  “  "€12*  "M12'  ~  <20b) 

Therefore,  for  the  dissipative  -  anisotropic  medium  the  orthogonality  relation  and 
dyadic  Green's  functions  can  be  expressed  in  the  following  way 

[i  ^  x  -o  ~  ^1-0  X  dS  =  2NJaP  <21) 


If  there  is  reflection  symmetry  in  the  waveguide  then 
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can  exist  independently.  This  approximation  is  also  valid  when  (ka)  <<1,  where 
2a  is  the  diameter  of  the  anisotropic  column. 

For  static  approximations  it  is  easy  to  show  that 


e  a'i 
.  a  _  *  1 

Y)  1  2? 

1  1  € 

r  r 


,  a 

** *  *  77 

r  r 


€*2 

z 


S  »  0 

M  a#  -  (t7t22  -  rj'i2) 
rj  sf  - 


(2) 


Now  using  these  relations  and  the  assumption  kjlft  <  <  1,  it  can  be  shown  from  the 
dispersion  relation  (1)  that  for  the  H-type  mode  the  right-hand  side  of  the  expres¬ 
sion  (1)  vanishes  and  for  the  E-type  mode  the  denominator  of  the  left-hand  side  of 
(1)  vanishes.  For  a  magnetic  current  ring  source  an  E-type  mode  will  be  excited 
in  this  static-limit  situation.  Whereas  for  an  electric  dipole  source  an  H-type  mode 
can  be  excited  in  the  static  limit.  It  is  of  practical  interest  to  consider  the  E-type 
modes  in  a  plasma  and  the  H-type  modes  in  a  ferrite.  The  following  dispersion 
relations  for  these  two  limiting  cases  can  be  expressed  in  the  following  way: 
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(n'ia)er 


J'(n'ia) 

n 

J  (rj'ia) 
n 


-  nef 


(p  (*a)K  (*b)  - 1  (ffb)K'  (/a) 
if  \  11  n  n  - 

*  -  wa;e2  |-j  ^b)R  ^a)  _  j  ^a)R 

Ln  n  n  n  j 

for  E-type  modes 


and 


J'  (u'ja) 

(an’j)wr  /gjjy-  W  • 


P  (ifc)K'  (tfa)  -  r  (tfa)K'  Kb) 

■n  n _ a  a  J 

I'(rfb)K  tta)-I  (*a)K'  (*b) 

„  n  n  n  n 

for  H-type  modes 


When  the  anisotropic  medium  completely  fills  the  waveguide,  i.  e. ,  when 
a  =  b,  the  above  two  relations  reduce  to  the  following: 


Jn(r}'ia)  =  0  ,  for  n  /  0 
and  for  n  =  0 ,  Jjiri\a)  =  0 


for  E-type  modes  (5) 


It  should  be  noted  also  that  in  this  case  a  change  in  sign  of  n  does  not  effect  any 
result. 


(n'ta)M 


J^(rj'aa) 
r  J  (rj’ja) 


=  om’ 


and 


for  n  =  0 ,  Ji(n'3a)  =  0 


f  for  H-type  modes  (6) 


j 


When  the  radius  of  the  waveguide  b  — ►  oo,  the  above  equations  (3)  and  (4) 
reduce  to  the  following  simple  forms: 

+These  results  agree  with  those  obtained  by  Trivelplece  ^8  J  except  for  a  change 
in  sign  in  the  term  nju'  of  equation  (4). 
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+For  a  ferrite  problem  only,  the  relation  n\  =  rj'j  cannot  be  satisfied. 
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For  n  =  0,  the  above  relation  becomes 

^rjVa^/^Jpj^ialJ^rjlal-arjilj^rjiaKJfCrjla^]  GQ(a)>/(a) 
e^[ezr>^(a)|2J1(T7'ia)JJnla)|^2-^r)i9(^(nia)+*I?(’liaJ^+2€2r7'i3Co(a)JJ(nla)) 

=  nMzGo(a){2Jo(nla)Jx(n\a)  (ni^sJ-arjlS^nlal+J^rjiaj)]  ^rj^SJa^rjla) 

(10) 

For  n  =  0,  the  dispersion  relations  (551)  and  (561)  can  also  be  derived  from  the 
relation  (1). 

It  should  be  mentioned  here  that  all  the  following  dispersion  relations  for 

various  special  cases  can  be  derived  from  any  of  the  general  relations  (551),  (561), 

and  (1)  of  this  Appendix.  The  purpose  of  writing  these  various  forms  of  the  general 

dispersion  relation  is  that  it  is  found  more  convenient  to  use  one  particular  form 

rather  than  another  for  some  special  cases. 

When  the  anisotropic  medium  completely  fills  the  w  *<%guide  (i.  e. ,  when 

a  =  b),  the  corresponding  dispersion  relation  (without  using  any  approximation) 

becomes  (since  G  (b)  =  C  (b)  =  ,  -  0 ,  (b)  =  0  =  S  (b),  at  a  =  b): 

n  n  or  rjhJ  n  n 

fto^l  fjlfo'ia)" 

*VrMl'iVi  [jn(„',a)J  -  tornfo 

(11) 
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r7*2S  Ji(rj'2a)Jo(r?'ia) 

"^M  =  J  (n’aa)Ji(n'ia)  ' 


for  n  =  0 


If  one  analyzes  equations  (341)  and  (351)  or  equation  (32a)  of  Appendix  A, 
it  is  easy  to  show  that  in  the  limit  e*  =  a*'  =  a^  — ►  0 


n'i*  = 


e  a'j 
z  1 


.  VL  .  V4 

^  ^  €  U  U 

r  r  r 


S  =  0 


M  =  -  (rj'a*  -  n'A  eza4  -  =  er£ 


-  =0,  c  a  -  c  rj'i2  =  e  M 

a^  z  4  r  * 1  r 


MrcVia  -  a'aez 


=  -*?  -f- 

r 


Now  if  one  divides  the  numerators  of  both  sides  of  (541)  by  and  uses 
relations  (13),  the  dispersion  relation  can  be  shown  to  have  the  following  form. 


after  rearranging  terms: 
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e  rm'iJ'  (ri'ia) 
r  1  n  1 

«2a4Cn(a  f 

MrnT1'2Jnn'2a^ 

vA***’ 

w«a4V 

Jn(n'ia) 

J  (iV2a) 
l.  n 

*  G  (a) 

n  .J 

karj 

For  an  isotropic  medium,  the  dispersion  relation  can  be  obtained  from  (14) 
letting  rj’i2  =  r?' 2  =  a4  =  ^€r~X2,  and  er  =  €z, 

For  axially  symmetric  fields  (i.e.,  when  the  ring  source  is  of  constant 

g 

strength,  "  _  •  0),  n  =  0,  and  one  obtains  from  equation  (14)  the  following 

ou 

dispersion  relations 


Jo(rj'ia) 

and 

jyjnVMn’* a) 

JQ(r)'2a) 


eaa4Co(a> 

y(a) 

o 


for  E-type  mode 


*<2a4s0<a> 

G  (a) 
o 


,  for  H-type  mode 


(15a) 


(15b) 


The  relations  (15a)  and  (15b)  can  also  be  obtained  from  (551)  with  appropriate 
limiting  procedures. 

It  should  be  noted  here  that  for  axially  symmetric  fields  and  e'  *  0  *  n\ 
E-type  modes  and  H-type  modes  can  exist  separately.  But  in  the  present  problem 
where  the  source  is  a  magnetic  current  ring  source,  only  E-type  modes  will  be 
excited  for  e'  =  0  =  n' .  Further  it  may  be  stated  that  even  in  an  isotropic  medium 
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g 

(i.  e. ,  c'  =  0  =  /u'  and  c  =  c  and  n  =  ju  )  if  ~rr  /  0,  E-type  and  H-type  modes 

rz  rz  dd 

cannot  exist  independently. 

If  n  /  0,  but  a  =  b  and  e'  =  0  *  n\  it  can  be  shown  either  from  (11)  or  from 
(14)  that  the  dispersion  relations  become 


Jn(n'ia)  =  0  , 

for  E-type  modes 

(16a) 

and 

J^n'aa)  =  0  , 

for  H-type  modes 

(16b) 

It  is  now  trivial  to  see  that  for  n  =  0,  a  =  b  and  e' 

=  0  =  /uf,  one  obtains  the 

following  dispersion  relations 

•Mrj'ia)  =  0  , 

for  E-type  modes 

(17a) 

Ji(n'aa)  =  0  • 

for  H-type  modes 

(17b) 

Relations  (16)  and  (17)  are  valid  for  isotropic  media  as  well  as  diagonally 
anisotropic  media.  However,  in  our  present  problem,  we  consider  only  (16a)  and 
(17a),  restricting  our  consideration  to  E-type  modes  (due  to  choice  of  the  source). 

+When  an  infinite  column  of  an  anisotropic  medium  is  placed  in  another 
unbounded  isotropic  medium,  the  corresponding  dispersion  relation  for  surface 

+It  should  be  noted,  however,  that  to  obtain  radiated  fields  in  the  present  situation, 
this  limiting  process  (namely  b  -+  a>)  is  not  valid. 
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waves  can  be  obtained  either  from  (531)  or  from  (541)  with  b  — ►  oo.  To  study  sur¬ 
face  wave  propagation  (which  is  also  slow  wave),  for  which  ^t/k  >  1,  one  can 
show  that  T)  is  purely  imaginary.  So  putting  rj  =  -j  &  ,  5  >0,  one  obtains 


J  (a)  =  -(2/*)fl  C$b)K  (5a)  -  I  (Sa)K  (5b)l 
n  L  n  n  n  n  J 

(18a) 

S  (a)  =  -(2/jr).fr(^a)K'(^b)  -  I'(5b)K»(5a)"| 
n  lq  q  n  n  J 

(18b) 

c  (a)  =  ( j2/jr)-fr(5a)K  (5b)  -  I  (5b)K»(5a)l 
n  l.  n  n  n  n  -J 

(18c) 

G  (a)  =  (j2/jr)  •  (~I'(5b)K  (5a)  -  I  <5a)K'(5b)l 
n  ^  n  n  n  n  J 

(18d) 

for  convenience  let 

J  (r)  =  I  (5b)K  (5r)  -  I  (5r)K  (5b) 
n  n  n  n  n 

(19a) 

S  (r)  =  I'  (5r)K*  (5b)  -  P (5b)K'  (5 r) 
n  an  n  a 

(19b) 

C  (r)  =  I'($r)K  (5b)  -I  ($b)K»(5r) 
n  n  a  an 

(19c) 

G  (r)  =  I*(5b)K  (5r)  -  I  (5r)K'(5b) 
n  n  n  an 

(19d) 

Also  it  is  easy  to  show  that  (if  n  <  $  b,  8  may  be  finite  or  very  large), 

J  K(5a)e*b 

V  — 7== -  ,  as  b  — f  co  (20a) 

“  J2*  $  b 
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When  n  =  0,  the  above  relation  becomes 

K  (iaJJ^n'ial+rj'i/iaKi^alJ  (n'ia)}  £e  K  (£a)Jx(f)'ia)+ear)'iKi(£a)J  (n't*) 

Ufa  K  (ialJ^n^al+n'aMjK^ialJ  <rj'2a)}  &c  K  (^aJJ^n'jaJ+ejn'jK^^aJJ  (rj'i*) 

C  Z  0  Os  z  o  o 

(22a) 

The  equation  (22a)  may  also  be  written  in  the  following  form: 


*  J1  (q,a*_T7,ia)  *3Lh,(ri'i2-‘n'i2) 

-  +  - J5 - Kja0a)J  (rj'ia)J  (rj'aa) 

O  0  o  o 


n'i  n'a 


(cjm  S  |  IcaM  S-*ij€  Ml 

-r,  Z  K  (*a)Kt(4a)J  (rtfa)  J^niaH  L  — z- ,  -J  k  (<5a)Kj(^a)J  (rj'2a)J1(r)\a)=0 
on  a  0  0  on  i  0  0 

(22b) 

When  the  anisotropic  medium  completely  fills  the  waveguide  in  such  a  way  that 
a  s  b  >)1,  the  corresponding  dispersion  relation*  can  be  obtained  from  (11)  letting 
a  >>  1.  It  should  be  noted  here  that  if  a  >>  1  in  equations  (21)  where  the  limit 
b  — *  oo  has  been  taken  the  result  will  be  different  from  that  obtained  from  (11)  with 
a  >>  1.  The  reason  is  that  if  a  >>  1  in  (11),  it  also  means  a  =  b  »  1,  but  if 
a  >>  1  in  (21),  it  means  that  b  »  1,  a  »  1,  yet  a  /  b.  The  difference  between 
the  above  two  results  can  be  shown  rather  easily  for  n  =  0  as  follows:  either  from 
(11)  or  (12)  one  obtains  the  following  dispersion  relation 


When  a  =  b  — >  oo,  there  is  no  dispersion  due  to  the  boundary. 
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r)'as  tan(i7'aa  -  ar/4) 
rj’jM  ~  tan(rj'ia  -  f/4) 


,  as  a  =  b  »  1 


and  from  (22a)  one  obtains 


s£f^Sin(r)ia-T /4)+rjiAJjCos(r]ia-T  /  4)}  ^e^Sin^a-T  /4)+e2r)\cos(rrta-x/  4) 

Sin(r)^a-x  / 4)+T7^cos(rj^a-ir7 4)]-  ”  ^e^Sinfrjia-*  / 4)+«jjrj^cos(r7ia-T  / 4) 

for  a  >>  1,  b  »  1,  but  a  ^  b.  (24) 


For  e1  *  0  *  /u1  =  a^,  the  dispersion  relation  (21)  becomes 


e  i  r)i J'(n\a)  f^K'^a)"]  Hu  4nfcJ' (t^a)  R' (jafl  ft 

r  n  ii  r  a  ^  ii  ^  ^ 

J  (n'ia)  K  (5  a)  J(rj'aa)  K  (5  a) 

n  J  ^  n  n  j  w 


n/(a+52) 

4 


Equation  (25)  can  also  be  obtained  from  (14)  directly  using  rj  =  -\b  ,  and  letting 


b — ►(». 


For  n  =  0,  the  E-type  and  H-type  modes  separate  and  one  obtains  the  fol¬ 
lowing  two  relations  from  (25) 


J/n'id 


eaa.Ki(ia) 

K  (5a) 
o 


for  E-type  modes  (26a) 


MjSnVMn'a*) 

J0(r)’2a) 


Ifea.K^a) 

K  (5  a) 
o 


for  H-type  modes  (26b) 


It  should  be  noted  here  that  for  isotropic  media,  r)\  and  rj'j  in  (25)  and  (26)  are 
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given  by 


r?’2  =  n'i2  =  n?2a  =  a4  =  ^Vr"**'  where  cr  58  cz  >r  =  \- 


For  e'  =  0  =  n\  n  =  0  and  also  a  =  b  »  1,  it  can  be  shown  either  from  (16)  or 


from  (23)  that  the  dispersion  relation  reduces  to 


cos(»j,1a  -  t/4)  =  0,  for  E-type  modes,  a  » 1 
sin(r)'aa  -  x/4)  =  0,  for  H- type  modes,  a»  1 


(27a) 


(27b) 


Now  it  is  also  natural  to  discuss  the  dispersion  relation  for  a«  1  and  b  finite. 
To  do  this  the  following  approximations  will  be  used  (see  jjL  7^J  ): 


J  (x)  —  ,  for  0<x  <<  1 

n  n  \ 

N  (x)  ~  (2/x)log  (fr~)  ,  for  0<  x  «  1 


where  or  =  1.781072 


Vx)  ~  "  •  n  *  0<x  «  1 

Kq(x)  -  -  log  (yJ,  0<  x  <<  1 
Kq(x)~  (|)  n  M,  0<  x  <<1 


(28a) 


(28b) 


(28c) 


(28d) 


(28e) 


J'  (x)  —  -  ,  n  >  1,  0  <  x  «  1 

n  2n(n-l) ! 
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N'(x) 

n 


nl  2 


jrx 


n+1 


,  for  n  ^  1,  0  <  x  1 


ni  2 


n-1 


K'(x) - X1 

n  n+1 


,  for  n^  1,  0  <  x  <<  1 


< 


It  will  be  assumed  that  r?a  « 1,  rj'ia  « 1,  rf2  a  « 1,  in  addition  to  a  <  <  1 . 
Then,  one  obtains  the  following  expressions: 


<a>~  N'(rjb)+  !  f— )  J'(rjb) 


n  «n  .  n  '  ir  \rfiL/  n 


2  n! 


J>>~-  isriii  T-^r)  J„w-  ^*7  Now 

'  2  n  ! 


Sn(a)  *W>  "  '  N'  (Tjb) 

“  ff/1  n  2n(n-l)!  n 


C  (a)~  J  (rjb)-  ^ 

“  jr(rja)  +1  n  2% 


,n-l 


2n(n- 1)  ! 


NW 


S  (a)  ~  - 
o 


N 

Nt(nb)  +  |  Ji(rjb)  log 
o(a)~  *  J0(r)b)  log  (-^)  -  NJnb) 

7^  Ji(nb)  -  “  ^(nb) 

7*  Jo‘*>  +  7  *.«■» 


for  n  >  1 


}(a)  ~  - 
(a) 


/s/ 


for  n  =  0 
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Thus,  for  a  «  1,  rja  «.  1,  a  <<•  1  and  a  <<  1,  various  forms  of  the 
dispersion  relation  can  be  obtained  merely  by  substituting  these  expressions  in  (541), 
(551)  etc. 
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APPENDIX  D 


A  NUMERICAL  EXAMPLE 


In  this  appendix,  results  of  numerical  calculations  will  be  presented  for 
a  special  case  of  the  physical  situation  described  in  (14a)  of  Chapter  n. 

A  normalized  value  of  the  electric  field,  E  ,  will  be  computed  as  a  function 

J* 

of  —  for  the  smallest  eigenvalue,  rj'i »  where  we  are  treating  the  case 


n’i2  >  0 
rft  <  0 


1  -  \  >  0 


1  -  y  >  0 


2  >  x  +  y  >  1 


with  the  specific  parameter  choices 


ka  ^  10 


e2  =  6  (which  represents  a  glass  tube) 
Hi  s  1 

tf/k  =  0  =  2.5  (>  JsH  ) 
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l 


a  6  =  ka  ^  =  ka 

,-3 


a  =  10 


c  =  b  =  2  x  10 


-3 


(2) 


First  the  lowest  eigenvalue  must  be  determined.  To  do  this,  first  some 
necessary  constants  must  be  computed  as  follows: 


-  c. 


G.  =  c  = 


I  (6b)  K  (6a)  -  I  (6a)  K  (6b)  >  0 
0  0  0  0 


It(6b)  K  (6a)  +  I  (6a)  iq(6 b)  ■?>  0 
o  o 


S  =  c_  =  ^(Sb)  Ki(6a)  -  li(6a)  iq(6b)  >  0 


c  =  I  Ua)  K  («b)  +  I  (6b)  Kt  (5a)  >  0 
4  i  O  O 


(3) 


r}\  (y)  will  then  be  computed  from  the  following  expression 


r*8 


-  ^2(1  -  y)  |^/(1  -  x)  +  y|  -  t(/xy  +  y  f(y)J,  (4a) 


2(x  +  y  -  1) 


where 


and 


^  =  IS2-!, 


^(y)  =  ip2  x2  +  4^  x(l  -  y). 


(4b) 


(4c) 


For  a  range  of  values  of  y  (where  1  -  x  <  y  •<  1). 

Starting  with  y  corresponding  to  the  smallest  value  of  r)\  obtained  from 
(4a),  the  function  G(rj'i ,  y)  given  below  will  be  computed  for  neighboring  y  until 
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a  change  in  sign  is  obtained.  £i.e.,  we  want  to  find  lowest  value  of  n'i/k  which 
satisfies  G(rj\  ,  y)  =  0 J 


G(rj’i  ,  y)  =  (p  +h'i2) 


(l-ylCjC^^all^a)  CjCgC^J^rjia)! 


n'i  P 


a)IQ(pa)  " 
2 


Jo(rj'ia)  Ij(pa) 

€P 


[e,  MCjOj  -  (1  -  y>sV3] 


r  -j  I  (pa)J1(rj,ia) 

L<*  -  -  e»sc4c2J  — 7^ —  *  0 


(5) 


where 


M 

S 


P2( y)  r)'aa  y  f(y)  rj\*(y) 

k2  k2  x+y-1  k2 

iil  -  y>d  -  +  2 

x  +  y-1  [  1  - x  PJ 

M-(p2+  T?'ia) 


(6a) 


(6b) 


It  should  be  noted  that  neither  r)\  =  0,  nor  p  =  0,  can  be  a  solution  to  (5). 

When  the  lowest  eigenvalue  i)\  is  determined,  |  |  will  be  calculated  as 

a  function  of  r  from  the  following  expression  where  all  relevant  parameters  are 
evaluated  at  the  previously  determined  lowest  tj'i  • 
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2<Je0e2  e  5,  2  cP  (r) 


‘2Z 

Jt= 1 


,  for  a  4  r  £  b 


F<  ^ 


where 


?i  =  ~ 


p  [Sc^t^a)  +  n'lCg^^'iaj) 
n'l  r4c2Ij(pa)  +  pc3IQ(pa)' 


»  [s€iIQ(pa)  -  MJjn'ta)] 

?a  =  71 - T-9 -  (8b) 

2^0)  +r}'ia) 

F1  *  L  er  €z  n'l*  P2  MR  +  umqp2 R'kjSt']  (8c) 

F2  =  "  J  L  er  e2  ^  P  ?1  [W  p2T’  k  &  c'  +  €rRn'ia  s]  (8d) 

F3  =  -iLVzp^  ?1  [w^R'k^'-e^Ml],  (8e) 

F4  =  '  L  Vz p2  ?  2  [u  V'1*  T'  k  06'  +  erT  rj'x*  s]  ,  (8f> 

P5  =  ^  [«„«*  «/{£«»  +  P  ^«b>}  +  (‘05*3{l?#b>+£  K?(*b))]  (8g) 

F6  =  ^  [v*V4b>{JIotfb)+;  K„“w}-i('o53I><*b>{il,£‘<‘W+Il«b>}] 


F  =  F 
7  6  * 
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p8  *  7^  Lv^  ^  (4b)  +  “o  *J  ’?<Sb)]  • 

jrue  c  e'kf) 

So  =  - I"  J  (n'xa)  -  ?x  I  (pa)  ]  , 

3  2erc2(pJV)  L  0  0  J 

AVe'k^ 

4|p2nla(p2+rjla)l2 


R»  = 


k3e  M(e  -  /33)  -  k^c  e'2 
r  r _ z 

we  c  cf 
o  z 

k  0  [erM  -  cgk2(€r  -  ft2)] 


k^  c  €'2  -  k2e  Sfe  -  02) 
z _ r  r 

UC  €  €? 
o  z 

*0  [erS  -  k2tz(tr  -03 


J(r)  =  I  (6b)  K  (6  r)  -  I  (Sr)  K  (6b) , 
o  00  0  0 


Jv.* 


€  =  -77—  (inM.K.S.  unit) 

O  OOf 


H  -  4j  x  10~  (in  M.  K.  S.  unit) 
o 
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ez  =  i-y  , 


e  = 


x  +  y  -  1 
x  -  1 


■I  S 


Ii  * 


y  >nr 

X—  1  * 


I.  = 


a 

2 


[J?(n’ia)  -  ^7  V^aKr^a)  +  J*o  (n'ia)j , 

0  Ja  —  fpJi^ia^a)  +  rj’i  Ii(pa)J,(n'ia)  1 , 

(p2  +  n;2)  L  J 

P2I  (pa)  Ii(pa) 

(pa)+  - - -  -  £<pa>J. 


■4-p 
-  v  p 


(6b)  + 


(6  a)  + 


2Mpa)  Ii(pa) 
pa 

2IQ(6b)  l!(6b) 
6b 

21  (6a)  Ufa) 
o _ 

6a 


(9j) 

(9k) 

(91) 

(10a) 

(10b) 

(10c) 


-  i  <*•>>] 

-  i  <«•>] 


(10d) 


i„  = 


6  7 


b2  r 


L(6b)K,(6b)  +  I  (6b)K  (6b)  + 
o  o 


1^(6  b)K1(6  bJ-1^6  b)KQ(»b) 
6b 


] 


2  r 

-  v[» 


I1(6a)K1(Sa)  +  Io(«a)Ko(6a)  + 


Io(«a)K1(6a)  -  ^(CajKjfa) 


6a 


o  ’  ~ 

I 


(lOef 
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*8  '  V 


2K  ((  b)Ki(5b) 


-  K  («b) 
o 


,  2K  (« a)K1(6a)  _ 

'(Sa)-  ~  °-  -  -K‘(8a) 

O  H.  O 


Discussion  of  the  Computation  Procedure  and  Results 


Starting  with  y  values  near  1.0,  it  is  found  that  n'i  increases  very  rapidly 

as  1  -  y  increases.  An  analysis  of  the  expression  for  G  shows  that  the  first  zero 

of  G  will  occur  shortly  after  J^rj'ja)  changes  sign,  i.  e. ,  when  r)\  a  is  slightly 

3  -3 

greater  than  2. 4,  i.  e. ,  n\  >  2.4  x  10  (with  a  -  10  ).  With  the  given  parameters 
it  is  found  y  is  very  near  1 . 0. 

After  finding  rj'i ,  and  the  corresponding  x  and  y,  for  which  G  =  0,  various 
expressions  in  equations  (8),  (9)  and  (10)  have  been  computed.  Since  the  parameters 
8  a,  6  b  are  large  in  this  calculation,  the  asymptotic  formulas  for  evaluating 
Iq(z),  Ii(z),  Kq(z)  and  K^z)  have  been  used  (where  z  stands  for  either  fa  or  8  b). 
With  the  parameter  used,  Kq(6  b)  and  Kj(fib)  are  very  small  and  nearly  all  terms 
involving  them  were  insignificant.  However,  the  form  used  included  all  operations 
and  can  be  used,  without  alteration,  for  any  set  of  parameters . 

It  should  be  noted  that,  although  some  of  the  Fg  and  If  in  (8)  and  (10) 
are  real,  some  pure  imaginary,  and  some  complex,  all  imaginary  terms  in  the 


summation  cancel  out,  and  hence 


Fj  1^  is  real.  This  fact  is  in  agree- 
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meat  with  the  theory  as  discussed  in  connection  with  the  equations  (82)  -  (84)  of 
Chapter  I,  where  it  is  stated  that  the  power-flow  in  a  non-dissipative  medium  is 
real. 

Computation  has  been  completed  for  the  following  cases: 

Casel)  x  =  0.7andx  =  0.5,  with  ka  =  2  x  10^  and  /3  =  2. 5(or  -  6. 25) 

Case  n)  x  =  0 . 7  and  x  =  0.5,  with  ka  =  10^,  f?  -  6.000025 

Casein)  x  =  0. 7  andx  =  0.5,  with  ka  =  7,  02  =  6.005102041 

The  following  pages  show  tables  of  values  of  many  of  the  variables  involved, 

c|E  I  c|EJ 

and  the  values  of  - — ,  for  0  <  a  <  b.  Also  graphs  of  - -  in  the  range 

m  m 

0  <  r/a  <  1.0  corresponding  to  the  above  cases  have  been  shown.  For  the  range 

«lE«l 

r  2  a,  the  values  of  — are  too  small  to  show  on  the  graph. 

The  behavior  of  all  the  graphs  plotted  here  is  more  or  less  the  same.  It  is 
the  nature  of  the  slow  waves.  The  higher  the  value  of  j3  =  ,  the  slower  the 

wave.  Moreover,  the  higher  the  value  of  e2,  there  is  a  minimum  value  of  f$,  for 
which  a  corresponding  slaw-surface  wave  can  propagate.  Since  in  the  above  com- 
uptation  e2  is  chosen  to  be  6,  the  minimum  value  of  0  is  greater  than  2.45.  On 
the  other  hand  a  larger  value  of  e2  will  permit  a  lesser  slow  wave  to  propagate. 
Although  the  above  statements  show  that  the  degree  of  slowness  of  the  surface 
waves  is  markedly  influenced  by  the  value  of  e2  and  hence  /3,  the  strength  or 
amplitude,  however,  of  these  waves  depends  on  various  other  parameters.  For 
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example,  in  the  Cases  I  and  m  the  graphs  show  that  higher  the  value  of  x  *  —  , 

O  E  | 

the  lower  the  amplitudes  of  — 1 —  .  On  the  other  hand,  in  the  Case  n  the 

m  CIE*I 

graphs  show  that  the  higher  the  value  of  x,  the  larger  the  amplitudes  of  - , 


although  the  values  of  0  in  all  of  the  above  cases  are  of  the  same  order.  Moreover, 

°  l^al  299 

the  graphs  of  the  Case  n  show  that  the  amplitudes  of  — 1 —  is  about  10 

m 

times  higher  than  that  of  the  Case  I  and  is  about  10  times  higher  than  that  of  the 

Case  m.  Therefore,  the  above  discussions  of  the  numerical  results  suggest  that 

for  any  practical  purposes  the  results  of  Case  n  will  be  of  greater  significance. 

All  the  graphs  plotted  here  change  monotomlcally,  because  of  the  higher 

value  of  0.  On  the  other  hand,  if  J3  is  small  (and  hence  the  smaller  value  of  the 

parameters  (pa) ,  6  a  etc. ),  it  is  expected  that  there  will  be  a  few  oscillations  of 
c  I  E_| 


— 1 - in  the  range  0  <  r  <  a.  In  support  of  this  statement  reference  may  be  made 

m  7 

elsewhere  [7}  . 

Finally,  it  should  be  noted  here  that  is  is  the  value  of  S  a  and  6  b  which 
played  the  significant  role  in  producing  tremendous  difference  of  amplitudes  of 

c  lE,l 

— 1 —  in  Case  I  and  either  Case  II  or  Case  in.  In  the  former  case,  the  value 
m 

of  6  is  much  higher  than  that  of  either  of  the  later  cases. 
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i 

t 


i 


i 


\ 


CASE  I 


Summary  of  Values  for  ak 


2xl02.  0  ‘  2.5 


Values  of 


X 

.7 

.5 

El 

For  x  * . 7 

For  x  = . 5 

y 

f 

*1 

p 

M 

S 

«2 

€r 

R 

T 

R* 

T» 

53 

L 

% 

a 

h 

h 

h 

k 

*5 

R*F6 
lm  F$ 

? 

% 

Vi 

*2*2 

*3*9 

r4*4 

f5% 

*•*6% 

Fgle 

Ef.Ib 

Coeff  1 
Coeff  2 

.9990731125 
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Summary  of  Value*  for  ak  ■  102,  &z 
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